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ABSTRACT 

We present Spitzer /IRS spectra of 31 TTS and IRAM/1.3mm observations for 34 
low- and intermediate-mass stars in the Cep 0B2 region. Including our previously pub- 
lished data, we analyze 56 TTS and the 3 intermediate-mass stars with silicate features 
in Tr 37 (~4 Myr) and NGC 7160 (~12 Myr). The silicate emission features are well re- 
produced with a mixture of amorphous (with olivine, forsterite, and silica stoichiometry) 
and crystalline grains (forsterite, enstatite). We explore grain size and disk structure 
using radiative transfer disk models, finding that most objects have suffered substantial 
evolution (grain growth, settling). About half of the disks show inside-out evolution, 
with either dust-cleared inner holes or a radially-dependent dust distribution, typically 
with larger grains and more settling in the innermost disk. The typical strong silicate 
features require nevertheless the presence of small dust grains, and could be explained 
by differential settling according to grain size, anomalous dust distributions, and/or 
optically thin dust populations within disk gaps. M-type stars tend to have weaker 
silicate emission and steeper SEDs than K-type objects. The inferred low dust masses 
are in a strong contrast with the relatively high gas accretion rates, suggesting global 
grain growth and/or an anomalous gas to dust ratio. Transition disks (TD) in the Cep 
0B2 region display strongly processed grains, suggesting that they are dominated by 
dust evolution and settling. Finally, the presence of rare but remarkable disks with 
strong accretion at old ages reveals that some very massive disks may still survive to 
grain growth, gravitational instabilities, and planet formation. 

Subject headings: stars: pre-main sequence - protoplanetary disks - planetary systems: 
formation - stars:late-type 
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1. Introduction 

Evolved protoplanetary disks are among the most interesting objects to understand how (and 
if) planet formation occurs in the disks around solar-type stars. Although there is well-accepted 
evidence suggesting that most IR excesses from protoplanetary disks have dissipated by ~6-10 Myr 
(e.g. Haisch ct al. 2001; Sicilia-Aguilar et al. 2006a, from now on SA06; Hernandez ct al. 2007), 

the way disks become optically thin and disperse after a (brief or long) transition phase is not 
clear. Different disk dispersal processes have been proposed, and observational evidence suggests 
that they all happen to some extent, but the relative importance and the interplay between them 
is unknown. 

Grain growth/settling to the midplane without fragmentation would render the disk optically 
thin in a timescale too short to be consistent with observations (Dullemond &; Dominik 2005; 
Brauer et al. 2008). Nevertheless, evidence for grain growth to millimeter-sized particles is found 
in (sub)millimeter observations as flatter dust emissivity curves and reduced disk masses compared 
to gas observations (Andrews & WiUiams 2005, 2007; Rodmann et al. 2006; Natta et al. 2007). 
High accretion rates observed in objects with negligible near-IR excess also point to strong grain 
growth/settling (Sicilia-Aguilar et al. 2010). Planet formation would substantially affect the disk 
structure, removing large amounts of mass and affecting the viscous evolution and accretion (Najita 
et al. 2007). UV/X-ray photoevaporation by the host star may not be able to fully disperse the 
disk in a reasonable time if it is too massive or too strongly accreting, but it can remove a large 
amount of mass, leading to inner holes and fast disk dispersal with the help of viscous evolution 
(Clarke et al. 2001; Alexander ct al. 2006; Gorti et al. 2009; Ercolano et al. 2009). In addition, the 
presence of stellar companions may contribute to the disk clearing (Bouwman et al. 2006; Ireland 
&; Kraus 2008), although it may not be the leading cause of inner holes (Pott et al. 2010). 

The observations of objects with no near-IR excess and strong mid-IR fluxes and thus optically 
thin/clean inner disks a few AU in size and optically thick outer disks, also called transition 
disks (TD; Skrutskie et al. 1990), suggests inside-out evolution to be an important disk dispersal 
mechanism, as predicted by Hayashi et al. (1985). Nevertheless, Spitzer revealed systems where 
the disk dispersal processes seem to operate over larger radial distances, resulting in disks with 
global dust depletion (Currie et al. 2009). The connection between the different disk dispersal 
mechanisms is complex. Photoevaporation is efficient at dispersing the full disk only once the disk 
mass and accretion rate have dropped below certain levels, and it is also favored by grain growth, 
which decreases the shielding of the gas to the FUV radiation (Gorti et al. 2009). Grain growth and 
settling are suppressed by gas turbulence (Schrapler &i Henning 2004), and may thus be important 
in the disk dead zones, but dead zones may disappear once the disk mass decreases, or never exist 
if the initial disk mass is too low (Hartmann ct al. 2006). In addition, the relative importance of 
the various dispersal processes is expected to vary with time (Alexander & Armitagc 2009). 

Although even the youngest star-forming regions contain evolved and transitional protoplan- 
etary disks (e.g. Hartmann et al. 2005; Fang et al. 2009), such objects with reduced IR excess 
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seem to become increasingly frequent with cluster age (Megeath et al. 2005; SA06; Lada et al. 
2006; Hernandez et al. 2007; Sicilia-Aguilar et al. 2009; Currie et al. 2009; Muzerolle et al. 2010). 
Nevertheless, there is still substantial debate regarding the definition of transition and evolved 
disks and where to place the limits to distinguish them from primordial ones, as quantifying mass 
depletion and the presence of inner holes is not easy without complete multiwavelength data and 
spatially resolved observations (Luhman et al. 2010; Currie & Sicilia-Aguilar 2011). The changes 
in the spectral energy distribution (SED) and thus, in the dusty disk component, observed with 
age are accompanied by a decrease in the accretion rate (Hartmann et al. 1998; Sicilia-Aguilar et 
al. 2006b, 2010), which even drops to undetectable or zero levels for roughly half of the transition 
objects (Sicilia-Aguilar et al. 2006b, 2010; Muzerolle et al. 2010). The optically thin emission of 
silicate grains (Calvet et al. 1992; Natta et al. 2000) also indicates dust evolution from ISM-like 
to larger and crystalline grains (Meeus et al. 2001; Bouwman et al. 2001; van Boekel et al. 2003; 
Kessler-Silacci et al. 2005; Furlan et al. 2005; Watson et al. 2009). Global disk structure (Lommen 
et al. 2010) and accretion/turbulence in the disk (Sicilia-Aguilar et al. 2007; from now on SA07) 
are thought to affect the grains we observe in the disk atmosphere (Dullemond & Dominik 2008; 
Zsom et al. 2011). Therefore, dust mineralogy can be used as a tracer of the energetic processes in 
the disk and the effect of turbulence and transport in the innermost regions. 

Evolved clusters and star-forming regions constitute an ideal laboratory to explore the char- 
acteristics and effects of disk dispersal, studying dust properties and disk structure for statistically 
significant samples of objects. One of the most interesting and well studied evolved regions is the 
Cep 0B2 association (Patel et al. 1995; Patel et al. 1998), located at 900 pc distance (Contreras 
et al. 2002) and containing the clusters Tr 37 and NGC 7160 (Platais et al. 1998). These clus- 
ters have been extensively studied at optical wavelengths in order to determine the membership, 
spectral types, extinction, the presence of accretion, and the accretion rate for a large number 
of solar-type members. The approximate median ages derived from the solar-type members are 
~4 and ~12 Myr for Tr 37 and NGC 7160, respectively (Siciha-Aguilar et al. 2004, 2005). Our 
Spitzer observations (SA06; SA07) revealed strong disk evolution. The disk fractions for the solar- 
type stars are ~45% in Tr 37 and ~4% in NGC 7160, and a large fraction of the objects present 
IRAC fluxes systematically lower than younger regions (e.g. Taurus; Hartmann et al. 2005) that 
suggest important grain growth/dust settling. Although among the objects with MIPS detections 
the number of evolved disks in Tr 37 is not much different than what is seen in Taurus (Luhman 
et al. 2010), the cluster lacks very massive and strongly accreting objects (with GM Cep being the 
only one in the region, Sicilia-Aguilar et al. 2008a) and a substantial fraction of objects have very 
low IRAC excesses and MIPS fluxes below the detection limits, characteristics of transitional or 
very evolved (strongly settled and likely dust depleted) disks. Up to 20% of the disks have no or 
very small near-IR excesses at A <6/im ([3.6]-[4.5]<0.2), being thus empirically classified as Tpj^. 



Note that there is still lack of consensus regarding the definition of TD. We apply here the definition we used 
previously in SA06, while some authors include as TD those showing global or homologous dust depletion (e.g. Currie 
et al. 2009), or use stronger near-IR constrains to ensure the presence of a cleared inner hole (e.g. Luhman et al. 
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Accretion rates or strong constraints on the accretion (via high-resolution Ha spectroscopy) are 
available for a large number of members (Sicilia-Aguilar et al. 2005, 2006b, 2010), and previous IRS 
observations revealed a large variety in the shape of the silicate emission and a weak correlation 
between grain size and accretion (SA07). 

Here we present 31 new IR spectra of T Tauri stars (TTS) in the Cep 0B2 region obtained 
with Spitzer IRS, and 1.3 millimeter observations obtained for 34 objects with the IRAM 30m 
telescope, both intended to study and constrain the structure of evolved disks. In Section [2l we 
introduce the new observations and also reanalyze our previous IRS/IRAC/MIPS Spitzer results. 
In Section [3l we model the properties of the silicate features and the global SEDs in order to derive 
the mineralogy and disk structure. We then discuss the different implications of the observed dust 
mineralogy and disk structure for disk dispersal evolution in Section [H summarizing the results in 
Section [5l 

2. Observations and data reduction 
2.1. IRS data 

We obtained 5.4 — 35 ^m low-resolution spectra {R = 60 — 120) of 31 members of the clusters 
Tr 37 and NGC 7160, within the Cep 0B2 association, with the Infrared Spectrograph (IRS, Houck 
et al. 2004) on-board of Spitzer. Our data were obtained within our Spitzer GO program 30523, 
observed between 2006 and 2008. To improve the statistics, adding up to 59 objects (56 of them 
solar-type stars), we also include here the data from our program 8 (already published; SA07), and 
five objects from program 3137, already public, (14-141, 72-875, 72-1427, 21364762, 213625070). 
The objects in program 30523 were selected among the more than 100 disked objects detected with 
IRAC and/or MIPS by SA06, avoiding those that presented a complicated background or nearby 
bright sources. We did not intend to obtain a complete collection, but to sample the different types 
of disk morphology (from flared, Taurus-like disks, to TD with inner holes) and accretion rates 
(ranging from ~10^'' to few times 10^^^ Mq/jv, and including some non-accreting TD) revealed 
by our Spitzer photometry and optical study. The complete sample with the basic properties of 
the objects is listed in Table [TJ Many of the disks are substantially evolved, with IR fluxes well 
below the Taurus median, although we also include a few very active disks, one of which is the 
variable star 13-277 (also known as GM Cep). CM Cep, together with the brightest among the TD, 
21392541, were also observed with the high resolution IRS modules. The targets for programs 8 and 
30523 were positioned in the slit using the high accuracy PCRS peak-up before the observations. 



2010). We also put our focus on discussing the different disk structures and the physical processes to which they can 
be related, rather than on the naming conventions. 

^The ID of objects listed in this work correspond to the numbering used in our previous optical papers (see 
Sicilia-Aguilar et al. 2005; SA06) 
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Our reduced low-resolution spectra are based on the droopres products processed through 
the S15.3.0 (program 30523) and S18.6.0 (program 3137) version of the Spitzer data pipeline. For 
consistency, we also re-reduced the data from program 8 using a more recent pipeline, S16.1.0, 
finding no significant differences with our previous data reductiorlfl. Pixels flagged by the data 
pipeline or by visual inspection of the images as "bad" were replaced with a value interpolated from 
an 8 pixel perimeter surrounding the errant pixel. Our data were further processed using spectral 
extraction tools developed for the "Formation and Evolution of Planetary Systems" (FEPS) Spitzer 
science legacy project (see also Bouwman et al. 2008), based on the SMART software package (Higdon 
et al. 2004). The spectra were extracted using a 6.0 pixel and 5.0 pixel fixed-width aperture in the 
spatial dimension for the observations with the first order of the short- (5.4 — 14 /um) and the long- 
wavelength (14 — 35 ;um) modules, respectively. The low-level fringing at wavelengths > 20 fim was 
removed using the irsf inge package (Lahuis & Boogert 2003). To remove any effect of pointing 
offsets, we matched orders based on the point spread function of the IRS instrument, correcting 
for possible fiux losses (see Swain et al. 2008 for further details). This was especially important 
for the objects from program 3137, which systematically presented offsets up to 2" from the real 
position. The background was subtracted using associated pairs of imaged spectra from the two 
nodded positions along the slit for most of the objects, also eliminating stray light contamination 
and anomalous dark currents. A few stars are embedded in nebulosity (73-472, 14-141, 21-998, 
21364762), having a background that suffers from variations on small spatial scales. For them, 
the standard method (subtracting the nodded spectra from each other) could not be used for the 
background correction, and we applied a customized routine for background subtraction. After bad 
pixel correction and flatfielding, the two-dimensional spectra for each of the nod-positions were 
extracted. The source positions in both nods were matched (they are typically at 1/3 and 2/3 of the 
spatial dimension) and the spectra combined. The source is then in the center of the spatial axis of 
the combined spectrum. For each wavelength resolution element a second- or third-oder polynomial 
was fitted to the regions to the left and right of the source spectrum, carefully avoiding the edges 
of the detector. These polynomials are evaluated in the central (source) region and subtracted 
as background. While this procedure results in increased uncertainty for the faintest sources, the 
extracted SEDs are in excellent agreement with the IRAC and MIPS photometry. 

The spectra were calibrated using a spectral response function derived from multiple IRS 
spectra of the calibration star t]^ Doradus and a marcs stellar model provided by the Spitzer 
Science Center. The spectra of the calibration target were extracted in an identical way as our 
science targets. The relative errors between spectral points within one order are dominated by the 
noise on each individual point and not by the calibration. We estimate a relative flux calibration 
across an order of ~ 1 % and an absolute calibration error between orders/modules of ~ 3 %, which 
is mainly due to uncertainties in the scaling of the marcs model. 

The high-resolution data for GM Cep and 21392541 were reduced using the standard pipeline 



Reduced spectra can be provided upon request by contacting the authors. 



- 6 - 



developed for the C2D data (Lahuis et al. 2006, 2007). The C2D pipehne is based on the SMART 
routines (Higdon et al. 2004), and includes advanced reduction and calibration routines developed 
for the C2D and FEPS data in order to perform the full aperture and optimal PSF extraction and 
pointing correction. The observations included separate sky observations, in order to improve the 
background subtraction. The data reduction starts with the SSC pipeline RSC products, which are 
corrected for crosstalk, in addition to the other corrections of the Basic Calibrated Data (BCD). The 
spectrum is extracted by using an optimal, wavelength-dependent PSF fitting. The full aperture 
spectrum is also extracted separately, in order to check for extended emission. The flux calibration 
is based on a large number of high-S/N calibrator stars and the MARCS stellar model provided by 
the Spitzer Science Center. The extracted spectra were corrected for fringing (Lahuis & Boogert 
2003), the orders were matched, and the data were corrected for the flux loses due to pointing 
offsets. The bright star CM Cep is well detected with the high-resolution IRS modules, but the 
object 21392541 is too faint, and its spectrum could not be extracted. 



2.2. IRAM 30m observations 

A total of 34 members of Cep 0B2 were observed at 1.3mm with the bolometer at the IRAM 
30m telescope (see Table [2]). The observations included most of the objects with IRS spectra, 
and took place within 3 different programs carried out between 2006 and 2009 (142-06, 132-07, 
145-08) as part of the MAMBO Pool Observations, using the MAMB02 117 channel bolometer. 
In addition, CM Cep had been previously observed with IRAM and the MAMBO 1 37 channel 
bolometer (Sicilia-Aguilar et al. 2008). 

All observations were performed as standard ON-OFF scans, with each scan including 10 min 
on-source and 10 min off-source integrations, divided in 1 min integrations following the ON-OFF- 
OFF-ON sequence. For each source, we obtained between 1 and 10 scans, varying the position 
(by observing at different times) and distance (30-70") of the off-source pointing for repeated 
observations. The beam of the telescope at 1.3mm is 11". While for most of the objects this is 
enough to avoid contamination by nearby sources or nebular emission, the bolometer observations 
revealed extended emission in a couple of cases (see Table [2]) , which are therefore excluded from 
this study. Skydip observations were performed every 2-3 hours or more frequently, depending on 
the weather, in order to determine the sky opacity, r. The pointing and flux calibration were done 
using the nearby sources LK Ha 234 and NGC 7538, from the IRAM Pool catalog, and the planet 
Uranus. 

The data were reduced using the MOPSIC pipeline, specially developed by R. Zylka for the 
reduction of ON-OFF observations. The MOPSIC pipeline uses the RSD (Resampled Signal and 
then the Difference) to subtract the atmospheric emission from the data while minimizing the 
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sky noisqj, which is measured from the correlated signal of adjacent pixels. The total emission is 
calculated taking into account the r measured in the skydip observations and the elevation of the 
source, and calibrated according to the flux calibrators above mentioned. The final result is the 
weighted mean of all the scans. For each source, we started reducing all individual scans, in order to 
detect potential problems. In particular, observations taken during the winter 2007 (program 132- 
07) suffered from sporadic technical problems due to malfunctions in the preamplifier. The affected 
scans, together with other scans suffering from poor weather or spikes, were labeled and removed. 
Finally, the good scans were reduced together for each source. The final fluxes (or 3a upper limits 
in case of non-detection) are given in Table [51 together with the estimated disk mass based on the 
Beckwith et al. (1990) prescription, for an opacity Ky=2 cm^ g~^ at 1.3mm. Interestingly, we do 
not find any correlation between the estimated disk masses and the accretion rates, except for the 
fact that the most massive disk corresponds to the strongest accretor, GM Cep. 



The IRAC and MIPS data for the cluster members had been presented in SA06. Given the 
improvement of the pipeline and calibration since then, we have re-reduced the data and redone the 
photometry for the objects analyzed here. In addition to our IRAC/MIPS data (AORs 3959040 
and 4316416 for Tr 37, 3959296 and 4320000 for NGC 7160), we have also analyzed the archive 
data from AORs 6051840, which include a few of the objects to the west of Tr 37 that were out of 
our original Spitzer field. 



The data reduction was done with the Mosaicking and Point-source Extraction pipeline (MOPE] 
The mosaics were created from the BCD data, using standard parameters for each band and in- 
cluding the masks for detector artifacts and associated status masks. The process included pointing 
refinement and background matching to produce uniform mosaics. Short and long exposures were 
treated separately. AORs 3959040 and 3959296 contained 5 individual images per pointing in both 
the short (Is) and long (26.4s) exposures, while the total exposure time of AOR 6051840 was 4.8s. 
MIPS AORs 4316416 and 4320000 mapped the region with median exposure times 78s. Aperture 
photometry was done with the APEX multiframe pipeline within MOPEX for IRAC channels 1,2, 
and 3. The presence of substantial extended emission in IRAC channel 4 and MIPS 24^m in Tr 
37 precluded the use of APEX for the photometry at these wavelengths, which was then done 
using IRAF standard tasks within the phot and apphot packages. The apertures and sky annuli 
were 3 and 12-20 pixels for IRAC, and 13" and 20-32" for MIPS, which are standard values. The 
corresponding aperture corrections were taken from the IRAC and MIPS Handbook, being 1.112, 
1.113, 1.125, 1.218, and 1.167 for the 3.6, 4.5, 5.8, 8.0, and 24/xm bands, respectively. Given that 
long-exposure images suffer more often from scattered-light artifacts and extended nebular emis- 



2.3. IRAC/MIPS 




^See |http://ir am.fr/IRAMFR/ARN/dec05/node9.html 

''http://ssc.5'pjtzer. caltech.edu/dataanalysistools/tools/mopex/ 
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sion, we preferentially use the short-exposure photometry for objects detected in the short-exposure 
mosaics with high signal to noise ratio (S/N). Nevertheless, in most cases the differences for objects 
detected in both short and long exposure mosaics are negligible. 

The newly reduced data (see Table [3]) agree very well with our previous results, and the use 
of smaller apertures improves the observations of the objects with complicated background within 
the Tr 37 globule, whose IRAC and MIPS magnitudes are now in better agreement with the IRS 
spectra. The new MIPS data for NGC 7160 also shows a better agreement with the IRS spectrum 
of 01-580, due to a better background treatment. The IRAC photometry of two objects (21364762 
and 72-875) is strongly contaminated by nebular emission, being thus excluded from the analysis. 

3. Analysis 

3.1. SEDs and silicate features 

For the solar-type stars that conform most of our sample, the region from ~5.4/im to ~35 
/im traces disk material located at ~0.1 to ~30 AU (Calvet et al. 1992; Chiang & Goldreich 
1997), and 1.3mm observations are sensitive to the bulk of the disk mass (Beckwith 1990). In order 
to complete our picture of the objects, we include the available optical photometry (see Sicilia- 
Aguilar et al. 2004, 2005, 2010). The spectral types and extinctions are known for most of the 
members (Sicilia-Aguilar et al. 2005). In the few cases where this information is not available, 
we take the initial extinction to be the cluster average (Ay=1.67ib0.45 mag; Sicilia-Aguilar et al. 
2005) and obtain an approximate spectral type and a corrected extinction by fitting the SED to 
different marcs models (Gustafsson et al. 2008). The SED fitting also lead to better extinction 
estimates for a few objects (72-875 and 21362507 have slightly higher extinctions than estimated 
before. Ay ~2 and 3.5 mag, respectively). The offsets between optic8l/2MASS/ Spitz er data in a 
few objects are typically related to known strong variability (e.g. 82-272, CM Cep). The SEDs, 
compared to the photospheric emission given by the marcs models (for the low-mass stars) and 
to the stellar photospheres in Kenyon & Hartmann (1995; for the intermediate-mass stars) are 
displayed in Figures [T]fT2] (arranged according to their SED types, see Section [3^ . Tables H] and [5] 
offer a summary of the previously published optical and 2MASS JHK data. 

Figures [13] to [TTl show the details of the 10;^m silicate features. Most of the sources present the 
characteristic silicate emission found in TTS, due to the presence of a mixture of small (0.1-6//m), 
warm (150-450 K) silicate grains with amorphous and crystalline components (e.g. Natta et al. 
2000; Bouwman et al. 2001; Meeus et al. 2003). Although silicate emission is present in most 
disks, it is not evident in a few cases. In addition to the intermediate-mass stars CCDM+5734 and 
KUN-196, the lO/um silicate emission is not visible in 72-875 nor in 11-2322. The first object is a 
TD with very low accretion rate (the lack of IRAC data does not allow to check the presence of 
near-IR excess, but the shape of the IRS spectrum suggest photospheric flux down to ~5-6 /im), so 
the lack of silicate emission could be due to a real depletion of small grains in the innermost disk. 
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On the other hand, 11-2322 is a Ml star with IR excess at all wavelengths, so the lack of silicate 

emission could be due to strong grain growth and settling (in fact, due to the S/N limitations, 
we cannot exclude the presence of some few-micron grains in this object), as has been suggested 
for the Coronet cluster (Sicilia-Aguilar et al. 2008b). A few more disks have weak (or no) silicate 
emission, including some TD (92-393, 24-1796) and normal disks (11-1209, 11-2131, 21-998, 21-33). 

A few other objects show atypical silicate emission. One of the TD (24-515) shows only 
evidence of crystalline silicates, suggesting that most of the amorphous grains are in large (few- 
microns) aggregates. The spectroscopic binary 82-272 presents prominent emission at ~11 /xm, 
suggestive of a dominating forsterite dust component, as has been observed for RECX-5 in the rj 
Cha cluster (Bouwman et al. 2010), and 21380979 could be another forsterite source. A few other 
sources have strong enstatite emission at ~9.3 ^m (e.g. 13-236, 21395813, 23-570). The emission 
of 21-2006 and 21384350 is mostly consistent with amorphous silica (note that 21384350 is also 
partially affected by the strong emission of a bright nearby object), and 12-2519 may be another 
silica source. Nevertheless, due to the complex shape of the continuum, extracting the silicate 
emission from the spectrum is complicated in TD and in objects with small near-IR excesses whose 
SEDs show a change in the slope around 6-12um (what we call "kink" disks, see SA06). The 
features observed in 24-1796 and 12-2519 could be affected by the strong change in the SED slope 
occurring at these wavelengths, although the crystalline features in 12-2519 reveal that at least 
some small crystalline grains are present in this disk. 



3.2. Disk mineralogy and grain size derived from the silicate feature 

In order to study the disk mineralogy, we fit the optically thin silicate emission using the 
Two Layers Temperature Distribution (TLTD) method, developed by Juhasz et al. (2009), from 
now on J09. We exclude from the fitting procedure and subsequent mineralogy analysis the objects 
with low S/N (S/N<20), those with no evidence of silicate emission, and those suffering from strong 
contamination by nearby sources or nebular emission in the 10/im region (e.g. 21384350). Although 
the silicate emission at ~10//m and ~20//m traces only a small part of the dust content in the disk 
(small grains in the optically thin disk atmosphere with temperatures T~150-450 K; Calvet et al. 
1992; Natta et al. 2000), the presence or absence of crystalline grains and the grain size can reveal 
information about heating processes in the disk, irradiation, transport and turbulence (Meeus et 
al. 2001; Honda et al. 2003; DuUemond et al. 2006; Dullemond & Dominik 2004, 2008; Watson et 
al. 2009: Zsom ct al. 2011). 

The TLTD method (sec J09 for details) reproduces the silicate feature using a multicomponent 
continuum, including emission from the stellar photosphere, an inner rim, and the optically thick 
disk midplane. It considers a collection of four different dust species with sizes 0.1, 1.5, and 
6.0 /xm: amorphous silicates with olivine and pyroxene stoichiometry (Dorschner et al. 1995), 
forsterite (Sogawa et al. 2006), amorphous silica (Henning &; Mutschke 1997); and, in addition, 
enstatite grains (Jager et al. 1998) with sizes 0.1 and 1.5 /xm (since the emission of 6/xm enstatite 
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grains is indistinguishable from the continuum in cases of low S/N; see J09). The mass absorption 
coefficients are calculated from the material optical constants assuming that the crystalline grains 
can be assimilated to a Distribution of Hollow Spheres (Min et al. 2005), and considering standard 
Mie theory for spherical particles for the amorphous dust. The main difference with standard one- 
or two-temperature fitting methods is that the TLTD model considers that the emission originates 
in regions characterized by a continuous distribution of temperatures. Therefore, the flux can be 
written as (see J09 for further details): 

. ^ , ^ /^J^ a, mill 27]- 2-qa 

Fu=F^^,ont + 2^2^A,iAti,j / -^B^{T)T dT (1) 
where the flux in the continuum is: 



Vri?^ /""^r.min Ott 2-qr 

F.,cont = Do-^B,{n) + Di —B,{T)T—dT (2) 

Tr.max 



r.max 

T 



L m.min 



'2ilT 2 — qm 

+ D2 I -^B,{T)T—dT. (3) 



Here, i?* and T^, are the radius and effective temperature of the star. The temperatures of the 
disk atmosphere, rim, and midplane (Ta, Ty, and T^) are parameterized as power laws of the radius 
(with exponents qa, qr, and qm, respectively) and vary between a minimum and a maximum value. 
The exponents of the temperatures (qa, qr, qm) and the coefficients of each contribution {Dq, Di, 
D2, and Dij) are fitted to the data using a genetic optimization algorithm (PIKAIA, Charbonneau 
1995; see J09). The uncertainties in the IRS spectra are estimated by adding random Gaussian 
noise to the original spectrum at the noise level, and then repeating the fit 100 times. The final 
silicate composition, the total crystalline fraction, and the average grain sizes are obtained as the 
average of the whole set, and the errors are derived from the standard deviation in the positive and 
negative direction. Since the mineralogy at diff'erent radial distances is not necessarily the same 
and the S/N of our spectra is also strongly wavelength-dependent, we fitted separately the 7-14 and 
17-35/im ranges, as we have done before for similar sets of data (Sicilia-Aguilar et al. 2008b). In 
the long-wavelength region (17-35/im), the emission of large amorphous grains strongly resembles 
the continuum, so the crystalline fractions may have a tendency to being overestimated. 

The results are displayed in Tables [6ll9] and the fits to the silicate features are shown in Figures 
[T3]|18i As in our previous studies, we find that the TLTD model with the assumed grain com- 
positions and size distribution reproduces very well the observed features at both short and long 
wavelengths. The residuals do not suggest the presence of other additional components, and the 
main source of uncertainty in the models is the S/N of the spectra. The sensitivity to large grains 
depends strongly on the S/N of the spectrum: The emission from grains larger than ~6/im is not 
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very different from the continuum emission in this wavelength range, so separating the emission 
of large grains from the continuum is only possible in objects with high S/N, and the maximum 
grain size detectable in noisy spectra (S/N~20-50) may be underestimated. Therefore, in order to 
extend the study of the grain size to all spectra, we also measured the peak of the silicate feature 
by normalizing the spectra to the continuum obtained by fitting a straight line between the average 
values in the regions 7.2-7.4^m and 12.6-12.7/im (the peak values are listed in Tabled]). This is 
a similar procedure to what we used in SA07, and provides a good estimate of the grain size in 
objects with low S/N. The normalized peak is strongly correlated with the grain size (see Bouwman 
et al. 2008), as objects with larger grains show shallower and broader features. For our sample, the 
correlation between measured silicate peak and fitted grain size is very strong, with a Spearman 
rank coefficient -0.57 (range -0.77 — -0.26) and a probability p=9xl0~^ that the two quantities are 
uncorrelated, even though the continuum fit used to measure the normalized peak is not the same 
defined by the TLTD models. The disk mineralogy derived here is further analyzed in Sections 14.11 
and 14.21 where we compare the dust characteristics with the stellar and disk properties. 



3.3. Tracing the global disk structure with RADMC 

In order to explore in more detail not only the disk mineralogy, but also the global disk 
structure, we have performed detailed modeling using a full 2D radiative transfer code. Our aim 
with this analysis is not to produce accurate fits to the SEDs of the individual objects in the 
region, but to understand the types of disk structures that reproduce the observed SEDs. Detailed, 
full-disk models offer an advantage over simple silicate-feature fitting procedure, which is that the 
whole disk structure is taken into account, allowing to explore the global disk details together with 
the grain size distributions. This is particularly important in objects where the SED slope changes 
at the wavelengths where the silicate emission appears, for which defining the continuum level with 
a simple model can be uncertain and lead to variations in the strength and width of the extracted 
silicate feature and thus errors in the estimated grain size. Given the various morphologies of the 
disks in our sample, up to 20% of our objects could be affected by this continuum uncertainty (see 
Figures [T3lfT71 and comments in Table E]) . 

We explored what are the minimal changes that need to be performed to a normal CTTS 
disk structure in order to reproduce the observed SEDs. For this purpose, we have selected 10 
prototype disks with different SED morphologies that represent the main types of objects observed 
in the Cep 0B2 region (see Table [TOl for a detailed summary of the prototypes and their properties). 
The prototype disks include: a normal TTS disk with silicate emission (model 1, prototype 11- 
2037), a very massive and flared disk (model 2, prototype GM Cep), a dust-depleted CTTS (model 
3, prototype 13-236), a strongly depleted and settled disk (model 4, prototype 11-1209), a "kink" 
or pre-transitional disk (model 5, prototype 13-125C0), a depleted but strongly accreting disk with 



^^This object had been previously classified as TD, but the new IRAC photometry produces a color slightly over 
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strong silicate emission (model 6, prototype the old disk around 01-580 in NGC 7160), a classical 
(large 24//m flux or "turn-up") TD with weak silicate emission (model 7, prototype 24-1796), a 
classical TD with strong silicate emission (model 8, prototype 14-11), a TD with no silicate and 
low 24/Lim flux or settled TD (model 9, prototype 92-393), and a settled TD with strong silicate 
emission (model 10, prototype 13-52). 

All disks, with some variations due to difl'erences in stellar luminosity and mineralogy, can 
be classified within one of these subtypes. Note that we exclude from this classification the 3 
intermediate-mass stars (CCDM-l-5734, KUN-196, and DG-481) since they are most likely TD 
with very large inner holes and low dust masses, or debris disks. The SEDs in Figures [T]fTT] have 
been arranged according to this classification, and the SED type is also listed in Table [1] for each 
object. Figure [19] helps to visualize the SED type, which is based on the SED shape and the silicate 
feature. We do not attempt to define a unique set of indices for SED classification, but to organize 
the SEDs observed in Cep 0B2 in order to explore the underlying physics and disk structures. 
Given that the disk slope is far from constant in the IR, especially for the disks with inside-out 
evolution, we checked the SED slope, a, at different wavelengths, considering for this the 2MASS, 
IRAC, MIPS, and IRS data. From the IRS data, we derive the fluxes at 5.8, 8, 12, and 17/im by 
integrating the flux in a 0.8 fim interval around these wavelengths. For longer wavelengths, we also 
derived the flux by integrating within a 2//m interval centered at 23 and 33/im. Although both the 
23 and 33//m photometric points could be affected by forsterite emission, strong crystalline features 
are only present in very few sources, so general trends should remain unaffected. In general, for the 
plots and the analysis, we give priority to the IRS-derived slopes, using the IRAC and MIPS values 
whenever no IRS data was available. In any case, IRAC, MIPS, and IRS data are in excellent 
agreement. The SED slope values are listed in Table [TTJ For the SED classiflcation, we note that 
we consider as TD those with [3.6]- [4.5] <0. 2 and small or negligible near-IR excess, since also a 
few normal CTTS (e.g. 14-183) have [3.6]-[4.5]<0.2. We also list as TD some objects with [3.6]- 
[4.5]'-0.2-0.25 and strong turn-ups (72-1427, 73-472, 21392541, 21392570), since their SEDs are 
more similar to our TD prototypes than to other models (e.g., "kink" disks or depleted CTTS). 

We used the RADMC radiative transfer code packag (Dullemond &: Dominik 2004) for 
modeling 3-D axisymmetric circumstellar dust conflgurations. This package was described and flrst 
used in Dullemond & Dominik (2004) , and has been extensively used and tested over the last 7 years 
for continuum radiative transfer in protoplanetary disks and for computing the disk structure (see 
for instance Pascucci et al. 2004; Pinte et al. 2009). Although the code assumes axially symmetric 
density distribution (2D), photon packages are followed in 3D. The code uses a variant of the 
Monte Carlo method of Bjorkman & Wood (2001) to compute how the stellar photons penetrate 
the disk and determine the dust temperature and scattering source function everywhere in the disk. 



our TD cutofT, being rather an example of a pre-transitional (Espaillat et al. 2010) or "kink" disk, according to its 
IR colors. 



See http: / /www.mpia.de /homes /duUemon /radtrans /radmc / 1 
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Using volume ray-tracing the spectra and images can then be determined at any inchnation. For 
evolved disks, which can be very geometrically thin, special care is taken of the grid, ensuring that 
the vertical structure of the disk is always well resolved. RADMC can handle multiple dust size 
and/or composition populations, each having its own spatial distribution. This makes the code well 
suited for studying how dust growth/fragmentation/settling affect the SED of the disk. In case of 
different dust distributions (e.g., for the innermost and outer disk), the temperature of each one is 
treated separately, although for simplification all the grains contained in a collisional equilibrium 
size distribution are assumed to have the same temperature, independently of their size. 

We run different RADMC models in order to reproduce the full SED of the prototype disks, 
using the stellar parameters effective temperature, stellar radius, and stellar mass Tg//, R*, 
derived from the spectral types and optical photometry (see Table [TO]l . The disks are modeled in 
the simplest way, assuming a pressure scale height (Hp) that varies as a power-law with the radius, 
Hp/RocR^/'' . We observed that models with well-mixed dust and gas in hydrostatic equilibrium 
failed to reproduce the disks, so the scale height at the outer disk radius, Urdisk/^disk, had to be 
set for each case. This procedure is equivalent to including some degree of settling in the model. 
Although settling is likely to affect differently grains of different sizes, our simple model assumes 
that all the grains have the same scale height and are well mixed. We also assume a negligible 
accretion rate, which is true for most of the objects in Cep 0B2, given their typical M~10~^- 
10~^MQ/yr rates. The disk surface density is estimated from the total disk mass, assuming that 
it follows a power law of the radius with exponent -1. In all cases, we take the gas-to-dust ratio 
to be 100, and the total disk mass is varied to fit the millimeter emission and the shape of the 
mid-IR flux as seen in the IRS spectra. This method allows to better reproduce the dust content in 
evolved disks where we may expect strong grain growth, anomalous gas to dust ratios, and variable 
viscosity and accretion rate with the radial distance. We note that for most objects, the total disk 
masses inferred from the millimeter and 20-30//m data are below the values one would expect to 
account for the typical accretion rates of lO^^-~^M0/yr, which is suggestive of strong grain growth 
and/or anomalous gas-to-dust ratios. The outer disk radius is relatively unconstrained, as the main 
parameter is the total assumed disk mass. For each case, we explored different sizes for the inner 
hole by setting the temperature of the disk inner rim, from no hole or dust distribution starting 
at the dust destruction radius, T~1500 K, up to few-AU cleared holes, or inner rim temperature 
~150-200 K. We explored different silicate grain size distributions, assuming a collisional power- 
law distribution with an exponent -3.5 and a minimum and a maximum grain size, which vary 
between 0.1-10/im and 100-10000/um, respectively. A total of 25% of amorphous carbon dust was 
also included in all the models, assuming the same grain size distribution for carbon as we do for 
silicate grains. Since the inclination of the objects is unknown, we assume always a 45 degrees 
angle. 

We regulated the amount of crystalline silicates in the disk by assuming that the innermost 
disk is fully crystalline up to a given (variable) radius (vcrys), after which the amount of crystals 
decays as predicted by simple mixing models (see Pavlyuchenkov & DuUemond 2007) until the 
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grains are totally amorphous. Given that our models are to be understood as a representation 
of a large collection of objects that have different silicate features, we did not try to reproduce 
exactly the silicate profiles of the prototypes nor to derive their mineralogy, but to obtain roughly 
similarly strong or weak features. Therefore, we gave priority to studying the amorphous grain 
population, excluding crystalline grains in the cases where their presence is not important/evident. 
For simplicity, we also do not assume any radial dependency of the chemical composition of the 
dust, since our analysis of the short- and long- wavelength regimes of the IRS spectra does not 
provide conclusive evidence for it. In most models, dust populations with different grain sizes 
for the innermost and outer disk are required in order to fit the SED. Here we followed a similar 
approach as for the crystalline grains, assuming a given uniform composition up to a radius iCcomp, 
after which it is progressively replaced by the outer disk composition (Pavlyuchenkov &: Dullemond 
2007). Table [T2] lists the properties of the best-fit models, which are displayed in Figures [20]|2H 



3.3.1. Model 1: Normal CTTS 

A normal CTTS disk can be easily reproduced with a relatively massive disk (Mrfj^fc ~0.01 M*, 
assuming a gas to dust ratio 100) with grains varying between 0.1 and 10000 //m. Large grains are 
required in order to fit the millimeter data within a reasonable disk mass. Nevertheless, we find 
that even such a disk presents some degree of settling, as the assumption of hydrostatic equilibrium 
between gas and dust invariably produces a too large near-IR excess to fit our prototype ij^l. This 
group is the most numerous among the objects observed here, as well as those with IRAC and MIPS 
data from SA06. This suggests that even the most pristine-like disks in the Cep 0B2 region have 
already experienced substantial grain growth and settling, which is consistent with the findings of 
Manoj et al. (2011) in the Chamaeleon I region. 



3.3.2. Model 2: Massive disk 

Reproducing a very massive disk like that of GM Cep does not lack difficulties. The relatively 
low far-IR flux, compared to the millimeter data, is hard to reconcile unless other elements are 
added to the model. In Sicilia-Aguilar et al. (2008a) we tried an extra black body-like component, 
that could be interpreted as an increase of flaring at some 30-50 AU. With RADMC, we manage 
to reproduce the strong near- and mid-IR flux with a very massive (0.052 Mq, assuming a gas 
to dust ratio 100) and strongly flaring or relatively thick disk (Hrdisfc/Rdisfc=0.2 at the outer disk 
radius) with a small grain population in the innermost disk (only 0.1 ^m grains) that turns into 
a mixed population (0.1-100 /im) after 3 AU. Nevertheless, this model still fails to reproduce the 
very strong silicate feature, and its flux in the mid-IR is higher than observed. Assuming a smaller 
disk mass lowers the mid-IR flux, but this fails to produce the high millimeter flux. Given that 
its high-resolution IRS spectrum suggest that the object is extended (see Section l474l) . the most 
plausible solution is the presence of an envelope or surrounding nebulosity heated by the star, in 
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which case the disk mass could be lower. Two more disks in TV 37 appear very flared as well, with 
SED slopes similar to GM Cep, 21374275 and 93-720. Nevertheless, the lack of millimeter detection 
for 93-720 suggests a much lower disk mass, and we have no evidence of the presence of an envelope 
for any of them, although 21374275 is embedded in a small patch of nebulosity near CCDM-l-5734. 

3.3.3. Model 3: Depleted CTTS 

Several objects in Tr 37 show very strong near-IR excesses, comparable to typical CTTS, but 
comparatively low mid-IR fluxes. Their silicate features are usually small (but noticeable). Such 
SEDs require a slightly flattened/settled disk, with vertical scale height Hrdisfc/Rdjsfc=0-13 and a 
normal dust population including submicron to millimeter grains, but the disk mass must be always 
extremely low (<2 10"'^ M,,,, assuming a gas to dust ratio 100). Such a low disk mass would be 
consistent with a classification within the homologously depleted TD (Currie et al. 2009; Currie 
& Sicilia-Aguilar 2011), although their near-IR excesses are strong and their accretion rates are in 
general high (~5x lO~^M0/yr for the prototype 13-236). This last point suggests that the total 
gas mass is probably much higher than what we estimate from dust observations. Some extra small 
grains in the innermost disk, not following the collisional distribution (see option B in Table [T2]l . 
a puffed inner rim, or differential settling, may account for the higher near-IR excesses observed, 
compared to a pure settled model, although too much differential settling would result in a silicate 
feature much stronger than observed (Dullemond & Dominik 2008; see also our model #6). A 
disk composed of only submicron grains, where grain coagulation has not been able to proceed, 
would also produce a similar SED. The very low mass of the disk and strong small dust depletion 
could also be indicative of planetesimal formation. In this case, the total gas mass of the disk 
could be higher, thus explaining the relatively large accretion rates observed. Note that we include 
within this group the star 21364762, given its steep IRS spectrum with very weak silicate features, 
although its flux at >14/xm is uncertain due to nebular emission, and from the high accretion rate 
suggested by its veiled spectrum (Sicilia-Aguilar et al. 2005), it could be that its disk is more 
massive than the rest of type #3 objects. 

3.3.4- Model 4: Globally settled/ depleted disk 

In order to reproduce this type of SED, with nearly zero near-IR excess and strongly settled 
at all wavelengths, a low disk mass (~0.005M*, assuming a gas to dust ratio 100), very small 
vertical pressure scale height (Jlrdisk/^disk=^-^'^)-, and relatively large minimum grain size (>6 
/im) are required. As in model #3, the low disk masses and settled SEDs make them candidates 
for homologous dust depletion (Currie et al. 2009), despite their full gas masses are hard to quantify. 
Nevertheless, the models that reproduce the steep slope and small excess in the near- and mid-IR 
usually generate too little flux at around ~5-6 lira.. Adding an extra population with small grains in 
the inner disk (as we do in models #5 and #10, see below) does not increase the flux in this region 
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due to self-shadowing, and increasing the disk thickness or flaring is not possible if we want to 
reproduce the mid-IR SED. Therefore, we speculate that these SEDs, strongly aff'ectcd by settling, 
may still have some variations in the vertical scale height with radius, a pufi"cd-up inner rim or gap 
edge, and/or some degree of differential settling resulting in different vertical scales depending on 
particle size. 

3.3.5. Model 5: "Kink" disk 

The presence of disks with low near-IR excess and a sudden change of slope around 6-12 /xm is 
characteristic of Tr 37, as appears reflected in the median disk SED (S A06) . The global structure 
of the disk can be well reproduced by a model with two different grain distributions, one for the 
innermost disk (<5 AU) with a large minimum grain size (~10 /im), and another one for the 
outer disk with smaller grains. A small vertical scale height is also necessary. Nevertheless, such 
a distribution cannot explain the strong silicate features generally observed. An extra small dust 
component in the innermost disk is needed. Instead, assuming that the innermost disk (<0.5-2 
AU) has a small grain population (e.g. 0.1 /xm grains) and the rest is populated by a coUisional 
distribution with larger grains (minimum grain size 1-10 //m) can reproduce the near-IR and silicate 
feature quite well (with a tendency to overestimating the near-IR flux), but fails to reproduce the 
high mid-IR excess and the "kink" shape. We therefore believe that these objects suffer from 
stronger grain growth in the inner disk than in the outer disk, but the (inner) dust distribution 
may be anomalous. They can be considered then as a special case of Model #9 (see below) with an 
extra small dust population in the innermost disk. An extra population of small grains has been 
also included to explain the SEDs of pre-transitional disks with inner gaps (Espaillat et al. 2010), 
which have similarly low near-IR excesses, strong silicate features, and strong excess at >6/xm. 

3.3.6. Model 6: Old surviving, dust-depleted disk with strong silicate 

This model is based on the strongly accreting (M~4x lO-^Mo/yr) but dust-depleted disk 
present around 01-580, the only accreting star in the 12 Myr-old cluster NGC 7160. The star 91- 
506 in Tr 37 has similar IR colors, although it is younger and its higher 24/xm flux and millimeter 
data indicate that it has a more massive disk. The SED is characterized by a low mid-IR excess as 
in the models of depleted disks, together with a strong near-IR and very strong silicate feature. A 
very low disk mass ('^5xlO~^M*, assuming a gas to dust ratio 100), together with strong settling 
(vertical pressure scale height Hr(ijsfc/Rdj6fc=0.09 ) are needed in order to produce the flattened SED, 
as in models #3 and #4. But, in contrast to the mentioned models, this prototype has a very strong 
silicate emission, which requires the presence of small grains (~0.1/nm) in a non- negligible portion 
of the inner disk (~2-3 AU). Alternatively, differential settling with a disk atmosphere populated 
only by small grains could also explain the observed SED. Indeed, DuUemond k. Dominik (2008) 
suggested that strong silicate emission in disks with weak mid-IR excesses would be a sign of strong 



-17- 



size-dependent dust sedimentation. As in cases #3 and #4, our model does not provide enough 

near-IR excess, which also calls for flaring or vertical scale changes between the innermost and outer 
disk, differential (size-dependent) settling, and/or presence of extra puffed-up walls around gaps. 
The lack of dust mass compared to the accretion rate may also hint the formation of planetesimals 
or planets within the disk. 

3.3.1. Model 7: Classical TD with weak silicate emission 

By classical TD we understand objects with reduced near-IR excess and 24/xm fluxes com- 
parable or higher than those of normal CTTS (Hernandez et al. 2007). Simple models with two 
populations of grains for the inner and outer disk cannot reproduce both the shape and silicate 
emission observed for our model ^7 (prototype 24-1796). In order to reproduce the general shape 
of the SED, we need a low disk mass and a low vertical scale height, in addition to a change in grain 
population around ~5 AU. Assuming a collisional distribution of grains l-lOOOO/nm in the inner- 
most 5 AU and a submicron (0.1 jim) dust population at longer distances, we produce a relatively 
good flt to the data, but the excess at 5-12 jiva is too high. A model with a physical dust hole (see 
model #8) also produces a good flt, although in this case the excess at 5-12 ^um is too low. We thus 
speculate that this SED type may result from a disk with either an optically thin dust population 
in its innermost part, or a change in flaring/thickness between the innermost and outermost disk. 
A secondary wall at few AU may also help to reproduce the large flux in the mid-IR, which can 
otherwise only be attained with a submicron dust grain population in the outer disk. 

3.3.8. Model 8: Classical TD with strong silicate emission 

In order to reproduce this type of disk, with prototype object 14-11, a physical hole with a 
minimum size 2 AU (tj„ <250-150 K) is required. Settled dust distributions with small grains 
are unable to reproduce both the negligible near-IR excess and the high mid-IR flux with a strong 
turn-up around ~8 iim.. The best results to reproduce the global shape of the spectrum are attained 
with holes ~3.5-6 AU in size, which are roughly in agreement with our toy model presented in SA07 
(consisting of a single black body with T~70 K to model an inner disk rim). In order to produce 
the observed silicate emission, an optically thin population of small grains needs to be included in 
the inner rim (or inside the hole) of the disk, as it has been done by Espaillat et al. (2010) for 
pre-transitional disks. A change in composition (from submicron grains to large grains) around 6 
AU produces good results in a model with a ~5.5 AU hole. The vertical scale height of the disk 
is harder to determine with the available data: the large mid-IR excess can be obtained using a 
relatively large vertical scale height {B.rdtsk/^disk^^-'^) with a large (~10 //m) minimum grain size, 
or a smaller ^rdisk/^disk with a dust distribution containing small grains (~0.1 /xm minimum grain 
size) . 
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3.3.9. Model 9: Settled TD with weak silicate emission 

The small near-IR excess and relatively flat SED of our prototype #9 (star 92-393) can be 
reproduced with a two-population grain model, including large grains in the innermost disk (20- 
10000 /im below 15 AU) and small grains only in the outermost part (0.1/im). Including a collisional 
distribution with small grains in this part fails to reproduce the high mid-IR flux, and a single, 
large-grain distribution also underestimates the mid-IR excess. No clean hole is required in this 
case, which is also in agreement with our observations of accretion in the SED prototype 92-393. 
The mid-IR bump could also be explained with a change in flaring/extra wall at 10-15 AU. This type 
of SED is thus a good example of inside-out evolution. The relatively low disk mass (5x10-3M^, 
assuming a gas to dust ratio of 100) could also be explained with planetesimal formation. The 
formation of planetesimals or planets in the innermost 10-15 AU would reduce the dust contain 
we see in this region, and could also cause the presence of a secondary wall or change in disk 
flaring/thickness at this distance. 

3.3.10. Model 10: Settled TD with strong silicate emission 

In order to reproduce the very small near-IR excess of a disk similar to our prototype #10 
(object 13-52) we need a very small vertical pressure scale height (HrdjsA:/R'disA:<0.03) in combina- 
tion with large grains (minimum grain size >1-10 /xm). Nevertheless, to produce the strong silicate 
feature we need to include some optically thin small dust (0.1/im). This small amount of dust is 
also enough to produce the tiny excess over the photosphere observed. Nevertheless, assuming a 
collisional distribution with grains as small as 0.1;um would produce too much near-IR excess and 
dilute the silicate feature, so the small grain component needs to be decoupled from the general 
dust distribution. A physical inner hole is not required in this case, as long as the main dust distri- 
bution has a large enough minimum size and the disk is very settled, although the lack of accretion 
in our prototype requires either a very low gas content or a physical hole or gap within the inner 
large-grain disk. This, together with the mass depletion inferred from the very low mid-IR fluxes, 
makes these type of disks good candidates to look for giant planet formation, where accretion onto 
a planet could explain the lack of significant M. 

4. Discussion: Structure and properties of disks in evolved clusters 

4.1. Dust evolution and stellar properties 

In this section, we contrast the results of our mineralogy analysis (Section 13. 2p with the known 
properties of the stars in our sample. Considering the limitations due to S/N and unresolved disk 
structure, we examine the correlation between different stellar and disk properties and the disk 
mineralogy inferred from the silicate feature. Several trends have been observed in the past. In our 
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previous study (SA07), we found a correlation between the strength of the sihcate peak (related 
to the grain size) and the age, in the sense that younger stars appeared to have weaker silicate 
peaks and thus larger grains in their disks atmospheres. This could be interpreted as an effect of 
turbulence stirring larger grains into the disk atmospheres of the younger and more turbulent stars, 
with older and more quiescent disks being more settled (such an effect has also been observed by 
Fang et al. 2011). We also found that the presence of flattened SEDs without silicate emission 
appeared more frequently among M-type stars than for K-type objects. 

Examining the new data, we confirm again the trend of the sources with later spectral types (M) 
to have weaker silicate emission among the solar- type stars (Figure [25]). Excluding the intermediate- 
mass (spectral types A and B) objects, which strongly differ in properties from the rest of the 
sample, we obtain a Spearman rank correlation coefficient of -0.37 (range -0.61 — -0.06), with 
a probability p=0.006 that the two quantities are uncorrelated. Given the correlation between 
silicate peak and grain size, this relation could be related to differences in disk structure, transport 
and/or grain processing depending on stellar mass, even though our sample contains objects down 
to spectral type M2.5 only. Weak silicate features have also been observed for later-type stars in 
different regions (Sicilia-Aguilar et al. 2008b, 2009; Furlan et al. 2011). The fact that in disks 
around colder, less luminous stars, a given wavelength traces a closer-in part of the disk than in 
more luminous objects could be partially responsible for these differences, but we emphasize that 
our sample is dominated by late-K and early-M stars, which are not so different in temperature 
and luminosity. In addition, we observe a tendency of M-type stars to have steeper SEDs, which 
has also been observed even in very young regions (e.g. the Coronet cluster; Sicilia-Aguilar et al. 
2008b, see also Muzerolle et al. 2010; Sziics et al. 2010). This could be caused by a larger degree 
of dust settling or grain growth in the innermost disk, although further evidence is needed from 
detailed modeling of a larger number of low-mass objects. In addition, the differences between 
solar-type and M-type stars are not evident in all regions (Sicilia-Aguilar et al. 2009; Furlan et al. 
2011), and other parameters (like the accretion rates) may also affect the SED appearance (Manoj 
et al. 2011). 

Although there is still a weak trend for weaker silicate features being more frequent among 
the youngest disks, its significance is weaker after adding the new sample (Spearman correlation 
coefficient 0.21, in the range -0.11 — 0.49, with probability p=0.14 that there is no correlation 
between both quantities; note that we exclude for this analysis the stars with spectral types G and 
earlier, which have uncertain ages). Our extended sample also lacks strong correlations between 
the normalized silicate peak or grain size and the accretion rate, which is related to the turbulence 
level (using the most recent accretion rate estimates from Sicilia-Aguilar et al. 2010). Performing 
a Spearman rank correlation test, the probability that the accretion rate and silicate peak/grain 
size are uncorrelated are 0.46 and 0.67, respectively. 

We had previously attributed these trends of silicate strength with age and accretion rate to the 
effect of turbulence, which is thought to be stronger at younger ages. (SA07). The main difference 
is that, instead of being focused on normal CTTS disks as that of SA07, our new dataset includes a 
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large number of objects with settled and evolved disks, including many disks with signs of inside-out 
evolution. Exploring the different RADMC disk models, we find that settling combined with grain 
growth may also increase the large grain population in the disk atmosphere, as it happens in cases 
with strong turbulence, so the net effect of turbulence may not be evident nor easy to quantify. 
Recent theoretical (Zsom et al. 2011) and observational (Fang et al. 2011) studies suggest that the 
connection between disk turbulence and the grain size in the disk atmosphere is real. 

When studying the age effects, we must also keep in mind that old disks (>8 Myr) are typically 
rare, so studying the disk mineralogy over the whole disk life span is always limited by the paucity 
of old objects. Moreover, given that half of the disks have dispersed by the age of ~4 Myr (SA06; 
Hernandez et al. 2007), the long-surviving disks are probably those with initially large masses 
(that can stand long periods of viscous evolution and photoevaporation) and/or those for which 
coagulation and/or settling are far less efficient than in the typical objects. In this last case, very 
old objects may not be representative of the typical processes occurring in most aging disks. 

Finally, as it has been observed before, we do not find any evident correlation between the 
crystalline fraction and the stellar properties. The fact that crystallinization and amorphization 
may occur at variable rates over the lifetime of the star (Abraham et al. 2009; Glauser et al. 2009) 
is probably the main cause of this lack of correlation. Since the example of EX-Lupi (Abraham 
et al. 2009) is the most remarkable one regarding strong and fast variation of grain properties, 
we have examined the mineralogy of the sources with known accretion variations (of 1 order of 
magnitude or more; Sicilia-Aguilar et al. 2010). Looking at these sources (72-875, 11-2146, 14- 
141, 11-2031, 13-1238, 82-272, GM Cep, 13-236, and 13-1250), we observe a tendency to having 
large (few microns) amorphous grains, but no evidence of any trend to higher or lower crystallinity 
fractions. We also do not observe any trend relating grain size and crystallinity fraction and, as 
had been previously found, crystalline grains are always small. 

4.2. Dust evolution and disk structure 

Here, we study the results of our mineralogy analysis (Section 13. 2p . contrasting them to global 
disk properties as derived from the SED slope at different wavelengths. Several previous studies 
have suggested that the silicate emission can be used as a tracer of the global disk structure. 
Bouwman et al. (2008) found a correlation between the amorphous silicate feature and the SED 
shape, suggestive of dust settling. Lommen et al. (2010) compared the millimeter emission with the 
strength and characteristics of the silicate feature, finding a weak correlation that could indicate 
global dust evolution throughout the disk, as has been suggested by Currie et al. (2009). Our IRS 
sample was especially selected to contain disks with very different structures, in order to study 
whether the grains in the disk atmosphere are tracers of more global processes affecting the whole 
structure of the disk. Many of the SEDs of low-mass disks in Cep 0B2 are substantially different 
from those seen in younger regions, showing evidence of dust settling (at a global scale or with 
a radial dependence, as in "kink" disks), global grain growth and dust depletion, and inside-out 
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evolution (from differential growth/settling to inner holes in the TD). Transition and "kink" disks 
(similar to the so-called pre-transitional; Espaillat et al. 2010), having a very low and/or negligible 
near-IR excess and an otherwise normal mid-IR excess, are prime candidates to look for inside-out 
disk evolution (as predicted by Hayashi et al. 1985) and to search for evidence of the different 
mechanisms that can evacuate the innermost disk of small dust (coagulation, planet formation, 
photoevaporation, binarity). On the other hand, fully settled disks and those with very low dust 
masses are suggestive of global or homologous disk evolution (Currie et al. 2009) happening in 
both the innermost and the outer disk. 

We have explored the disk mineralogy derived from the silicate feature (in particular, the flux 
at the normalized silicate peak and the grain size) as a function of disk structure. Figure [26] displays 
the SED slopes (a, see Section 13. 3p at different wavelengths together with the peak height of the 
normalized silicate feature. Note that, for most transitional and "kink" disks, the slope changes 
sign between 8 and 24 fim, due to the turn-up in the SED. We find that objects with the most 
settled inner disks (including TD) have low silicate peaks, which are typical of large grain sizes. A 
Spearman rank correlation test confirms a strong correlation between the peak height and the SED 
slope at any wavelength (see Table [T3|) . with probabilities that the two quantities are uncorrelated 
varying between 6xl0~^ and 0.04. The correlation is stronger for slopes including short and long 
wavelengths, which trace the disk structure over several AU. This points out that strong grain 
growth in the innermost disk is probably the main cause of the SED geometry of transitional and 
"kink" disks, and it also suggests that grain growth strongly affects the disk structure (for instance, 
due to increased settling). In addition, the accreting TD tend to have have comparatively weaker 
silicate peaks than the non-accreting ones. If their reduced IR excess is caused by grain growth 
in the innermost disk it could explain why those objects are still accreting, as grain growth does 
not need to affect the gas content (see disk model #9 and Figure [28]) . The presence of large (>1 
^m) and processed grains (with normal crystalline fractions) in all the TD, including those with 
SEDs consistent with inner holes, offers a strong contrast to the findings of Watson et al. (2009) 
in Taurus, where the silicate features of TD with inner holes are nearly pristine. This is again a 
confirmation that the TD that are found in clusters of different ages do not necessarily correspond 
to the same type of physical objects (Sicilia-Aguilar et al. 2010). If the presence of pristine grains 
in the Taurus TD is due to dust replenishment from the outer disk, it could also indicate changes in 
the transport mechanisms or in the global dust processing. Differences in the structure of TD with 
various ages is also to be expected taking into account the different timescales of the various disk 
clearing mechanisms (planet formation, grain growth, photoevaporation; Alexander &: Armitage 
2009). 

We have also explored the potential correlation between crystallinity and disk structure, finding 
no significant results (Table [T3|) . There is a marginal correlation (p=0.04-0.09 in the Spearman rank 
test) between the forsterite to enstatite ratio and the innermost disk structure (Figure [27l) . This 
is consistent with Bouwman et al. (2008), who suggested that young disks (which have typically 
large near-IR excesses) would tend to form only forsterite in their inner disks due to their non- 
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equilibrium conditions. The slopes at longer wavelengths do not show any correlation with the 
forsterite to enstatite ratio, so the effect seems to be exclusive of the innermost disk, although the 
poor S/N of most spectra at long wavelengths does not let us test this hypothesis. Nevertheless, 
the forsterite to enstatite ratio is never too large in Cep 0B2, so the trend should be tested on 
a larger sample. We do not find the weak trend of Watson et al. (2009) that there are more 
crystals at shorter wavelengths (innermost disk regions): For the 6 objects that have enough S/N 
at longer wavelengths, we find a similar crystallinity fraction, with small variations, although our 
sample is too small for statistically robust results. In addition, detecting crystalline grains at longer 
wavelengths is easier than tracing the amorphous ones, so our long-wavelength crystallinity fraction 
may be overestimated. 



4.3. Disk dispersal and survival and the nature of transition disks in Cep OB2 

From the analysis of the previous sections, in particular, the comparisons with the model SEDs 
obtained with RADMC (Section 13. 3p and the trends observed with mineralogy (Sections 14.11 and 
14. 2p . we can draw general results about disk evolution in intermediate-aged clusters. Although 
we do not have millimeter observations (or detections) for all disks, our models revealed that an 
important deal of information about disk mass (low or high) as well as dust distribution (presence 
of small or large grains) in the planet-formation area can be already inferred from the ~20-30^m 
region in the IRS spectra, as it had been previously found by Currie & Sicilia-Aguilar (2011). Of 
course, a better constrain on grain growth and grain distributions will be only provided by deeper 
millimeter /submillimeter observations, or by far-IR data from Herschel. 

We observe important evidence of dust and disk structure evolution. On one hand, the low 
millimeter fluxes and the presence of objects with low mid-IR excesses suggest that dust evolution 
and grain growth up to millimeter sizes (and probably larger) is generalized throughout the disk. 
Not only some disks require relatively large grains in the innermost part, but a substantial portion 
of the dust mass must be forming large grains, as expected from global dust evolution (Currie et 
al. 2009). The inferred disk masses are sometimes very small, assuming a normal gas to dust ratio 
(~100), even in cases where we still detect strong accretion (e.g. 11-2146, 13-236, 13-669, 82-272, 
12-1091, among others). We thus expect these disks to have far more gas than what their dust 
emission suggests, which would also explain the lack of correlation between the accretion rate and 
the dust mass derived from millimeter observations already mentioned in Section 12.21 This could 
indicate an anomalous gas to dust ratio, together with a probable large maximum grain size. In 
case of very dust depleted disks like 13-236, the very low dust masses required to fit the SEDs 
could even suggest the formation of planets or planetesimals within the disk. We also find that 
even the most pristine-like disks in Cep 0B2 require some settling in order to properly fit their 
SEDs (in fact, CM Cep is the only one that does not require settling but rather extra flaring). 
Such a generalized dust settling has also been observed even in younger regions, like Chamaeleon 
I (Manoj et al. 2011). 
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In addition to this global grain growth and settling, we also observe significant evidence of 
inside-out dust evolution. More than half of the disks in Tr 37 show strong variations in the SED 
slope between the innermost part (<6-10//m) and the rest, consistent with inside-out evolution 
(SA06). In some cases, the near-IR excess is low enough to empirically classify the object as transi- 
tional ([3.6]-[4.5]<0.2). Our modeling of selected objects proves that non-standard disk structures 
can reproduce the SEDs of "kink" and transitional disks. Objects with a strong change of slope 
("kink" SED) invariably require different grain populations in the innermost and outer disk, as 
long as the vertical scale height is assumed to run smoothly as a single power-law throughout the 
disk. Surprisingly, a few of these disks do not require the grains in the innermost disk to be larger: 
a non-collisional distribution populated by submicron grains alone in a settled innermost disk can 
reproduce the near-IR excess and silicate feature quite well. Nevertheless, the large "kinks" typi- 
cally require the inclusion of some several-micron grains. We can interpret this result if small and 
large grains are not well mixed, as could result if large grains are settled, leaving a disk atmosphere 
populated mostly by small dust. Other authors have also proposed adding a small dust population 
within a disk gap for such objects with very strong silicate features (Espaillat et al. 2010). Includ- 
ing such a small dust population could improve our fits in the cases where the general structure 
of the disk that reproduces the global SED shape fails to produce strong enough silicate features. 
Size-dependent dust settling (with most large grains being confined to the optically thick midplane 
of the disk) could also be invoked to increase the strength of the silicate emission (Dullemond & 
Dominik 2008). Conversely, inside-out grain growth and/or grain filtering at gap edges due to 
the presence of planets (Rice et al. 2006) could explain different grain populations between the 
inner and the outer disk. Grain filtering due to planet formation, resulting in small grains in the 
innermost disk, also has the interesting side effect of increasing the gas to dust ratio (Rice et al. 
2006), which could explain the relatively large accretion rates observed in Tr 37. 

Clean inner holes are required to reproduce the objects with smaller near-IR excesses and large 
turn-ups at longer wavelengths. Since many of these objects do not seem to be actively accreting (in 
contrast with the objects with near-IR excesses, all of which display indicators of ongoing accretion), 
the presence of a physical hole is indeed the best explanation. In fact, the strong silicate features 
observed in many of these non-accreting transition disks suggest that the disks still contain a large 
population of small grains, albeit at a distance larger than the dust destruction radius. The disks 
with inner holes and no accretion are of great relevance, since they could be produced by the recent 
formation of a (giant) planet at few AU (Lin & Papaloizou 1986, Bryden et al. 1999, Quillen et 
al 2004), which would disrupt the disk and prevent accretion by opening a gap. An alternative 
explanation for TD without accretion is photoevaporation by the central star (Alexander et al. 
2006), although the number of non-accreting TD in Tr 37 seems in conflict with the timescales 
predicted for photoevaporation (Gorti et al. 2009). On the other hand, strong grain growth in the 
innermost disk is required to explain the transitional disks with weak or no silicate features. This 
result agrees well with the fact that all of these objects still show active accretion (see Figure [28]) . so 
the innermost disk must still contain a large amount of gas. In both cases (with and without hole), 
the flux mid-IR region reveals either inside-out evolution (objects with high mid-IR flux require the 
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presence of small grains at a few AU distance) or global grain growth (in objects with very steep 
SEDs in the mid-IR). 

Finally, we find two puzzling objects in the region. One of them is the massive disk around 
GM Cep (see Section [4.41 for a detailed discussion). Although the high-resolution IRS observations 
suggest that there could be some remnant of an envelope, the high disk mass and the strong (and 
variable) accretion rate are anomalous for a star in a 4 Myr-old cluster. The other one is the 
star 01-580 in NGC 7160, which despite being located in a cluster with median age 12 Myr, has 
a very strong silicate feature, strong near-IR excess, and an accretion rate of several times 10^*^ 
M0/yr (Sicilia-Aguilar et al. 2005, 2010). Considering its accretion rate, the disk around 01- 
580 must have been more massive than the one around GM Cep and similar to that of GW Ori 
(Mathieu et al. 1995). Such objects pose strong constraints to the possibility of planet formation via 
gravitational instability (Boss 1997, 2003), and also suggest that strong and global grain growth may 
not always be efficient. Therefore, also these anomalous objects may provide important information 
to understand the physical processes involved in planet formation. 

4.4. GM Cep: Jets, photoevaporation, and extended emission 

Of the two objects observed with the high-resolution IRS modules, GM Cep and 21392541, 
only the former was detected. This object has the brightest mid-IR emission among all the TTS in 
Tr 37, being twice as bright as the Herbig Be star MVA-426 at 70 ^m (SA06). We examined the 
spectrum looking for gas lines and found evidence of [Ne II] 12.81;um, [Ne III] 15.56//m, and [Si II] 
34.82^m in emission about 5-7(7 over the continuum (see Figure [29l) in the PSF-extracted spectrum. 
Examining the full aperture spectrum, we also find marginal evidence of PAH emission at 11.3;um. 
We do not find evidence of molecular emission in the spectrum, but it must be noted that the strong 
mid-IR continuum emission of the source does not allow to exclude weak lines due to a contrast 
problem. We find significant differences between the PSF-extracted and full aperture spectra in all 
the [Ne II], [Ne III], and [Si II] lines (see Figure [29]). suggestive of extended emission. In addition to 
the differences in the extraction, [Ne III] and [Si II] are more typical of H II and photodissociation 
regions (Osterbrock 1989), which also points in the direction of nebular contamination near the 
source. The relatively low contrast of the lines against the strong continuum emission does not allow 
for further interpretation in terms of dynamics. So far, no extended material (cloud nor envelope) 
has been detected around GM Cep at millimeter (Patel et al. 1998) nor Spitzer wavelengths, 
but the presence of a remnant envelope or cloud material could explain the high millimeter fluxes 
observed for this source, as well as the difficulties in fitting its SED with a simple disk model 
(Sicilia-Aguilar et al. 2008a). Moreover, if a remnant envelope or cloud material are present, as it 
is seen around other relatively isolated objects in the cluster (e.g., around the group formed by the 
binary CCDM+5734 and the F9 star 21374275; SA06), it could indicate that GM Cep is younger 
than the main Tr 37 cluster, which would be more in agreement with its high and variable accretion 
rate. 
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[Ne II] emission requires X-ray/EUV ionization (Glassgold et al. 2007), which could be due 

to the normal stellar activity, to the presence of jets/winds, and/or to active photoevaporation 
(Pascucci et al. 2007; Lahuis et al. 2007; Hollenbach & Gorti 2009; Giidel et al. 2010; Pascucci 
&; Sterzik 2009). GM Cep was one of the Tr 37 stars detected in X-ray with Chandra (Mercer et 
al. 2009), with observed luminosity Lx=l.lxl0^° erg s~^ (3.25x10^° erg s~^ after correcting for 
absorption). The IRS spectrum reveals an integrated [Ne II] luminosity L[jve//]=2xl0^° erg s~^. 
Both hx and Lj^ve//] appear higher than the typical values found in CTTS (Giidel et al. 2010), 
except when sources with strong jets are considered. While no definite confirmation of an extended 
jet has been found for GM Cep, its optical spectra reveal the presence of a strong stellar wind 
(Sicilia-Aguilar et al. 2008a). The double-peaked O I emission at 8446A is also suggestive of a 
bipolar wind/outflow, so the strong [Ne II] emission observed in GM Cep could be related to jet 
emission, as it has been detected toward T Tau (van Boekel et al. 2009). Active photoevaporation 
and puffed-up disks have also been suggested for increased [Ne II] emission. The high mid-IR fluxes 
observed in GM Cep suggest a very fiared disk, but due to the strength of the [Ne II] emission 
(more than one order of magnitude of the value to be considered within the normal CTTS range 
for its Lx] Giidel et al. 2010), the presence of a jet is the most plausible candidate. 



5. Summary and conclusions 

We present new Spitzer /IRS spectra for 31 TTS and IRAM 1.3mm observations for 34 low- 
and intermediate-mass stars, all within the clusters Tr 37 and NGC 7160 in the Cep 0B2 region. 
The observations were selected to cover a representative sample of SEDs and disk types observed 
in Tr 37, together with the only 2 low-mass disks detected in NGC 7160 (SA06). Including our 
previously published IRS spectra (SA07), we analyze a total of 56 IRS spectra of solar-type stars 
plus 3 of intermediate- mass objects with silicate features. We study the new data, together with the 
previously published optical, near- and mid-IR observations, in order to examine the correlations 
between silicate dust features and the stellar properties, accretion, and disk structure. 

We analyze the silicate emission features with a TLTD method, finding that they are consistent 
with a mixture of amorphous silicates (with olivine and forsterite stoichiometry), amorphous silica, 
and crystalline forsterite and enstatite grains with sizes 0.1-6/im. No other components are required 
in order to reproduce the observed features, and in general, the amount of crystals in the disk does 
not exceed 20%, with a typical value being 5-10%. The lack of S/N at long wavelengths in most 
objects does not allow to explore the radial dependence of the crystallinity fraction nor chemical 
composition. 

A large number of disks present strong signs of dust evolution: settled SEDs, evidence of 
optically thin inner disks, presence of "kinks" or slope changes in the SEDs. We thus characterize 
the global disk structure by measuring the SED slope at different wavelengths, and create a series 
of prototype disks by means of the RADMC radiative transfer code. This allows us to determine 
what types of disk geometry and dust distributions reproduce the behaviors observed in our sample. 
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The comparison with the models reveal that even the most typical CTTS disks in Cep 0B2 
have already suffered strong grain growth and settling. This general requirement of strong settling 
has also been observed in younger regions (e.g. Chamacleon I; Manoj et al. 2011), and thus seems 
to occur very early in the disk lifetime. In addition, we find that about half of the disks present 
signs of inside-out evolution, seen as a variable SED slope at different wavelengths (with a large 
negative slope at short wavelengths, changing to a near-zero or positive slope at longer A). Such 
a SED geometry can be reproduced by radially dependent grain populations or by changes in the 
disk flaring/thickness with radius. 

We find that TD (especially, accreting ones) and "kink" disks appear strongly correlated to 
grain growth in the innermost disk and typically show flattened and small silicate features. Even 
the objects whose SEDS can be explained by inner holes show processed grains, in contrast with 
what is seen in the TD with holes in Taurus (Watson et al. 2009), which tend to have nearly pristine 
silicate features. This may indicate that TD are physically different objects at different ages, as 
has been previously suggested (Alexander & Armitage 2009; Sicilia-Aguilar et al. 2010). While TD 
in Taurus could be mostly related to giant planet formation (still subject to migration, due to the 
high mass of the disks), substantial grain growth and settling (and eventually, photoevaporation) 
can explain the older TD in Tr 37. It could also be a sign that, at older ages, the mechanisms 
replenishing the inner disk with pristine grains are not effective, or the outer disk has also suffered 
grain processing. 

In order to reproduce the near- and mid-IR SEDs of transitional and "kink" disks, very settled 
disk structures and, in some cases, radially-dependent dust distributions are required, in addition 
to typically very low disk masses {Maisk <5xl0""^ M*, considering a standard gas to dust ratio of 
100). This suggest that, although evolution seems to operate faster in the innermost disk, there 
is also substantial evolution and dust depletion at a more global scale, as it has been suggested 
by Currie et al. (2009) and Lommen et al. (2010). We also find that many of our SEDs require 
the presence of submicron-only (--^O.l/um) grains in the inner disk. This could be caused either 
by differential settling, with larger grains being preferentially in the midplane, or by a radially 
dependent grain population. Populating the inner disk with small dust grains could be seen as 
inconsistent with insidc-out evolution, but may be explained by dust filtering at the edges of the 
gaps formed by planets (Rice et al. 2006). 

The differences in disk structure in our varied sample dilute the observed correlation between 
accretion rate (or turbulence) and the silicate feature, although a correlation is still seen in samples 
of normal CTTS disks (SA07; Fang et al. 2011). The relatively strong accretion rates observed in 
the stars in the Cep 0B2 region (of the order of few times 10~^ M0/yr) are in contrast with the 
small dust masses inferred from millimeter data and radiative transfer models, which could be an 
indication of a larger-than-normal gas to (small) dust ratio. This also points out the importance 
of taking into account the gas content to determine the evolutionary status and the fate of disks, 
including the possibility of planet formation. 
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Finally, the IRS high-resolution spectrum obtained for the variable star GM Cep reveals that 
it is probably extended and thTis surrounded by cloud material or a remnant envelope, which could 
be responsible for its strong and variable accretion rate. Its IR [Ne II] emission is also consistent 
with the presence of jets, which is in agreement with the strong wind signatures observed in its 
optical spectrum. 
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Fig. 1. — SEDs for the objects with SED type ^^1, normal CTTS (see text). The optical, 2MASS, 
IRAC, and MIPS data are represented by magenta dots, with their errors (usually, smaller than the 
dots). The IRS spectra are shown in black with errors in cyan. Photospheric models (Gustafsson 
et al. 2008; Kenyon & Hartmann 1995) are shown in green. The data are corrected for extinction 
using the Ay derived for each object. 
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Fig. 2. — SEDs for the objects with sed type #1, normal CTTS (continued). Same as Figm'e[T] 
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Fig. 3. — SEDs for the objects with sed type #2, massive CTTS (see text). Same as Figured) 
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Fig. 4. — SEDs for the objects with sed type #3, depleted CTTS (see text). Same as Figured! 
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Fig. 5. — SEDs for the objects with sed type #4, globally settled/depleted CTTS (see text). Same 
as Figure [H 
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Fig. 6. — SEDs for the objects with sed type #5, "kink" disk (see text). Same as Figure [H 
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Fig. 7. — SEDs for the objects with sed type #6, dust-depleted disk with strong sihcate feature 
(see text). Same as Figured) 
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Fig. 8. — SEDs for the objects with sed type #7, classical TD with weak silicate (see text). Same 
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Fig. 9. — SEDs for the objects with sed type #8, classical TD with strong silicate (see text). Same 
as Figure m 
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Fig. 10. — SEDs for the objects with sed type ^^9, settled TD with weak sihcate (see text). Same 
as Figured! 
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Fig. 11. — SEDs for the objects with sed type #10, settled TD with strong sihcate (see text). Same 
as Figure [H 
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Fig. 12. — SEDs for the intermediate-mass objects, which are most hkely dust depleted TD with 
large holes or debris disks. Same as Figure [TJ 
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Fig. 13. — Silicate features at in the ~10^m region for the objects in Tr 37. For each object, the 
upper panel shows the IRS spectrum in black, with errors in green. The fitted TLTD model is 
shown (if appropriate) as a magenta line, and the corresponding continuum is depicted as a dashed 
line. The lower panel shows the residua (black) compared to the errors (green). 




Fig. 14. — Silicate features at in the ~10;um region for the objects in Tr 37 (continued). Same as 
Figure [m 
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Fig. 15. — Silicate features at in the ~10;um region for the objects in Tr 37 (continued). Same as 
Figure [ni 




Fig. 16. — Silicate features at in the ~10/xm region for the objects in Tr 37 (continued). Same as 
Figure [IS 
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Fig. 17. — Silicate features at in the ~10/im region for the objects in Tr 37 (continued) and NGC 
7160 (lower right panels including the objects 01-1152, 01-580, and DG-481). Same as Figure [T3l 
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Fig. 18. — Silicate features at in the ~20-30/im region for the objects in Tr 37 with good S/N 
at long wavelengths. The S/N in these spectra ahow to identify features (11-2046, CCDM-l-5734, 
82-272, 13-277, 93-720) as well as to exclude the presence of them (21395813, KUN-196). 
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Fig. 19. — SED slope and normalized silicate peak for visualization of the SED classification scheme 
(see text). The SED types are determined considering the whole SED and the relative strength 
of the 10;um silicate feature, here we present the SED slope at H-[5.8] vs. H-[24] (left) and the 
normalized silicate peak versus the inner disk slope at [3.6]- [5.8]. The correspondence of SED class 
and symbol is as follows: red circles = ^ 1, green circles = # 2, blue circles = ^^3, cyan circles 
= #4, magenta circles = ^7^5, yellow circles = #6, black circles = 7,^7, red triangles = 7^8, green 
triangles = #9, blue triangles = #10. 
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Fig. 20. — Models for normal (left, #1) and massive(right, #2) CTTS, compared to the SEDs and 
IRS spectra of their prototypes 11-2037 and GM Cep, respectively. 
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Fig. 21.— Models for depleted CTTS (right, #3) and globally settled/depleted disk (left, #4), 
compared to the SEDs and IRS spectra of their prototypes 13-236 and 11-1209, respectively. An 
upper limit at 1.3mm is shown as a down-pointing arrow. For Model 3, we display the results of 
assuming a uniform dust distribution (3A, dashed line) versus the case with small grains in the 
innermost disk (3B, bold line; see Table [T2] for details). 
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Fig. 22. — Models for "kink" disk (right, ^^5) and long-surviving CTTS disk (left, 7^6), compared 
to the SEDs and IRS spectra of their prototypes 13-1250 and 01-580, respectively. An upper limit 
at 1.3mm is shown as a down-pointing arrow. For Model 5, we display the two fits: A, with 
small grains inside, which reproduces well the silicate emission (bold line), and B, with large grains 
inside, which fits the general structure and the slope change of the disk (dashed line; see Table [12] 
for details about models 5 A and 5B). 
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Fig. 23. — Models for transition disks with strong mid-IR excesses, for disks with weak sihcate 
(left, #7) and strong silicate (right, #8) features, compared to the SEDs and IRS spectra of their 
prototypes 24-1796 and 14-11, respectively. For Model 8, we display the models with large grains 
and large vertical scale height (A, bold line) and the model with small grains and small vertical 
scale height (B, dashed line; see Table [T2] for more details about the models). 




Fig. 24. — Models for settled transition disks with no turn up and no silicate feature (left, #9) 
versus silicate feature (right, #10), compared to the SEDs and IRS spectra of their prototypes 
92-393 and 13-52, respectively.An upper limit at 1.3mm is shown as a down-pointing arrow. 
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Fig. 25. — Normalized silicate peak versus spectral type. The intermediate mass stars (KUN-196, 
CCDM+5734, DG-481, and also the F9 star 21374275) are marked with boxes and shifted by 20 

subtypes for a better display in the figure, with their spectral types being B9, A3+A5, A7, and F9, 
respectively. Stars with spectral types MO or later have normalized peak fluxes lower than earlier 
objects, as it had been suggested in SA07. The errors in the spectral types (not shown) are 0.5 to 
1 subtype. 
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Fig. 26. — Disk slope (a) at different wavelengths versus normalized silicate peak, which is a tracer 
of grain size. The different a trace the disk structure at different disk radii. There is a tendency 
of more settled and transitional disks to have lower silicate peaks (or absence of silicate emission) 
and thus larger grain sizes. This suggest the importance of grain growth for the innermost disk 
structure. The correspondence of SED class and symbol is as follows: red circles = # 1, green 
circles = # 2, blue circles = #3, cyan circles = #4, magenta circles = #5, yellow circles = #6, 
black circles = ^7, red triangles = ^8, green triangles = #9, blue triangles = #10. Intermediate- 
mass (spectral types B and A) objects are marked with a box. Among the vertical lines, the dashed 
line represents the slope of a bare photosphere (spectral type K7), the dotted line is the slope of a 
geometrically thin, optically thick disk (Kenyon & Hartmann 1987), and the dotted-dashed line is 
the slope for the median of the Taurus disk (Fur Ian et al. 2006). 
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Fig. 27. — Disk slope (a) at different wavelengths versus forsterite to enstatite ratio. The different 
a trace the disk structure at different disk radii. There is a weak tendency of very flared disks to 
have higher forsterite to enstatite ratios. The correspondence of SED class and symbol is as follows: 
red circles = # 1, green circles = # 2, blue circles = #3, cyan circles = #4, magenta circles = #5, 
yellow circles = #6, black circles = #7, red triangles = #8, green triangles = #9, blue triangles 
= #10. Intermediate-mass (spectral types B and A) objects are marked with a box. Among the 
vertical lines, the dashed line represents the slope of a bare photosphere (spectral type K7), the 
dotted line is the slope of a geometrically thin, optically thick disk (Kenyon & Hartmann 1987), 
and the dotted-dashed line is the slope for the median of the Taurus disk (Furlan et al. 2006). 
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Fig. 28. — Comparison of the SEDs (including silicate features) of an accreting and a non-accreting 
transition disk with similar spectral types. 
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Fig. 29. — Lines observed in the high-resolution IRS spectrum of 13-277 (GM Cep). The PSF- 
extracted spectrum is shown as a thick magenta histogram; the full aperture spectrum is displayed 
as a thin black histogram. The differences in both extractions could be indicative of extended 
emission. 



Table 1. Summary of IRS observations and stellar properties. 



J.CLL gCU 


rinnrrlinatp^f' 79000^1 


Mode 


Tnt Timpl'^^ 


Sp.T. 










Sil Ppak 


V /' ^llXXXlCXXLiO 


Tr 37 






















73-758 


21;35;08.39 -1-57:36:03.2 


Low 


59-59-19-19 


K6.5 


1.9 


1.8 


8 





1.8+0.1 




72-1427 


21:35:16.32 -1-57:28:23.0 


Low 


44-183 


Ml.O 


1.4 


3.3 


8 




1.7+0.2 




73-472 


21:35:18.61 -1-57:34:09.2 


Low 


1451-1451-975-975 


K5.0 


1.7 


3.3 


7 


A 


1.2+0.1 




73-311 


21:35:24.51 -1-57:33:01.1 


Low 


366-366-315-315 


M1.5 


1.1 


2.8 


4 


A 


1.5+0.1 




73-71 


21:35:30.21 -|-57:31:16.4 


Low 


1934-1934-975-975 


K6.0 


2.1 


1.8 


4 


2.1^ 


1.4+0.1 




72-875 


21:35:49.80 +57:24:04.6 


Low 


44-183-94-44 


M0.5 


2.0 


2.6 


7 


o.utUl/0 


1.0+0.2 


SBl:, variable M 


21362507 


21:36:25.08 +57:27:50.3 


Low 


44-183-94-44 


MO.O 


3.5 


26 


1 


0.09l«:J| 


1.9+0.1 


SBl 


21364762 


21:36:47.63 +57:29:54.1 


Low 


44-183-94-44 


K6.0 


3.1: 


5.6 


3 


A 


1.0+0.2 


Extended? 


14-141 


21:36:49.41 +57:31:22.0 


Low 


44-183 


K6.0 


1.7 


0.8 


3 


<2.3 


1.2+0.1 




11-2146 


21:36:57.69 +57:27:33.3 


Low 


13-13-13-13 


K6.0 


2.6 


2.5 


1 


0.99101 


2.3+0.1 




11-1209 


21:36:58.55 +57:23:26.1 


Low 


29-29-126-126 


K6.0 


0.6 


1.4 


4 


6+0* 
"■"-0.2 


1.0+0.2 




11-2322 


21:37:01.96 +57:28:22.3 


Low 


13-13-13-13 


Ml.O 


0.9 


0.7 


3 


1 6+28 


1.2+0.3 




11-2037 


21:37:07.07 +57:27:01.1 


Low 


44-44-13-13 


K4.5 


1.6 


2.5 


1 


<9 


2.0+0.1 


1 

ai 
00 


14-11 


21:37:10.30 +57:30:19.1 


Low 


122-122-13-13 


Ml. 5 


1.8 


1.4 


8 





1.2+0.3 


14-125 


21:37:10.54 +57:31:12.4 


Low 


366-366-189-189 


K5.0 


1.7 


3.0 


1 


o.3t°;? 


1.5+0.1 


1 


11-2131 


21:37:12.19 +57:27:26.6 


Low 


13-13-13-13 


K6.5 


2.3 


3.0 


3 


o.3tg;i 


1.1+0.2 




11-2031 


21:37:15.95 +57:26:59.5 


Low 


13-13-13-13 


K2.0 


1.7 


4.5 


1 


o.3tg| 


2.8+0.1 




14-160 


21:37:27.32 +57:31:29.5 


Low 


366-366-189-189 


K5.0 


1.8 


3.1 


1 


A 


1.6+0.1 




14-183 


21:37:38.53 +57:31:41.7 


Low 


44-44-13-13 


K6.0 


2.0 


4.8 


5 


16+0-12 

"■^"-0.07 


2.1+0.1 


SBl: 


CCDM-l-5734 


21:37:40.93 +57:33:37.5 


Low 


13-13-13-13 


A3+A5I 


I.7I 


4.0* 


D 




1.0+0.2 


Visual binary 


21374275 


21:37:42.75 +57:33:25.0 


Low 


117-88-88-88 


F9.0 


3.0 


4.0* 


2 


A 


1.6+0.1 




12-2519 


21:37:51.07 +57:27:50.2 


Low 


366-366-189-189 


K5.5 


1.6 


2.8 


4 




1.4+0.1 




12-1091 


21:37:57.62 +57:22:47.6 


Low 


366-366-189-189 


K2.0 


2.8 


27 


1 


<0.4 


1.8+0.1 




13-1238 


21:37:59.26 +57:36:16.2 


Low 


366-366-189-189 


Ml.O 


2.6 


1.8 


1 


6+1-1 


1.4+0.1 




82-272 


21:38:03.53 +57:41:35.4 


Low 


13-13-13-13 


G9.0 


3.6 


13 


1 


<0.9 


2.1+0.1 


SB2 


13-669 


21:38:09.28 +57:33:26.2 


Low 


488-366-189-189 


Kl.O 


2.2 


6.3 


1 


6+0-2 
"."-0.2 


2.2+0.1 




21380979 


21:38:09.79 +57:29:42.8 


Low 


1451-1451-731-731 


Ml.'' 


1.7/ 


4.0 


10 


PA 


1.9+0.1 




54-1781 


21:38:16.12 +57:19:35.7 


Low 


975-731-377-377 


Ml.O 


1.2 


2.5 


10 


Q + 0.5 


1.4+0.1 




13-1877 


21:38:17.03 +57:39:26.5 


Low 


366-366-252-252 


K7.0 


2.0 


2.3 


1 


1 6+1-2 


1.5+0.1 




13-277 


21:38:17.31 +57:31:22.0 


Low 


50-38-50-50 


G8.0 


2.5 


1.9 


2 




2.5+0.1 


Variable M 




21:38:17.31 +57:31:22.0 


High 


189-366 
















13-236 


21:38:27.42 +57:31:08.1 


Low 


366-366-189-189 


K2.0 


1.8 


5.4 


3 


51+0.24 
'J.3J--0.13 


1.2+0.1 




12-2113 


21:38:27.47 +57:27:21.3 


Low 


13-13-13-13 


K6.0 


2.5 


1.1 


1 


6*- 


2.4+0.1 


SBl 


13-157 


21:38:28.12 +57:30:47.1 


Low 


13-13-13-13 


K5.5 


1.2 


2.4 


1 


4 2+2-9 


1.9+0.1 




13-52 


21:38:32.55 +57:30:16.1 


Low 


1451-1451-731-731 


K7.0 


1.3 


4.1 


10 





1.7+0.1 




91-155 


21:38:34.75 +57:41:28.0 


Low 


44-44-13-13 


M2.5 


1.2 


1.7 


1 


PA 


1.6+0.1 





Table 1 — Continued 



Target 


Coordinates(J2000) 


Mode 


Int.Time(s) 


Sp.T. 


Av (mag) 


Age (Myr) 


SED type 


M (10 M0/yrj 


Sil. Peak 


21384350 


21 


:38:43.50 


-1-57:27:27.0 


Low 


1451-1451-731-731 


M2.0 


1.7 


7.0 


9 





1.4+0.2 


54-1547 


21 


:38:44.48 


-1-57:18:09.6 


Low 


59_59_19_19 


K5.5 


1.1 


9.9 


5 


0.18+°" 

—0.06 


2.7+0.1 


91-506 


21 


:38:58.07 


-1-57:43:34.3 


Low 


1451-1451-731-731 


K6.5 


1.4 


6.5 


6 


0.09+°°I 

—0.04 


2.8+0.1 


12-1617 


21 


:39:04.68 


+57:25:12.8 


Low 


1451-1451-731-731 


Ml.O 


1.6 


1.6 


5 


0.2+Of 
^ — 0.1 


1.7+0.1 


13-1250 


21 


:39:12.17 


+57:36:16.9 


Low 


44-44-13-13 


K4.5 


1.4 


3.9 


5 


3+°-] 

0.1 


3.0+0.1 


21392541 


21 


:39:25.41 


+57:33:20.2 


Low 


488-488-189-189 


K6-f 


1.7/ 


1.1 


8 


A 


2.9+0.1 




21 


:39:25.41 


+57:33:20.2 


High 


5780-7254 














21392570 


21 


:39:25.70 


+57:29:45.5 


Low 


1451-1451-189-189 


K6^ 


Ijf 


4.0 


7 


PA 


1.6+0.1 


21393104 


21 


:39:31.04 


+57:47:14.0 


Low 


366-366-189-189 


K2f 




2.6 


1 




3.5+0.1 


24-515 


21 


:39:34.10 


+57:33:32.1 


Low 


426-426-144-144 


M0.5 


1.1 


5.6 


8 


0.10+°i« 

—0.05 


1.1+0.1 


21-998 


21 


:39:34.80 


+57:23:27.7 


Low 


1451-1451-731-731 


K5.5 


1.9 


9.7 


3 


0.3+°? 


1.4+0.1 


21-33 


21 


:39;35.64 


+57:18:22.3 




366-366-489-489 


MO.O 


1.7 


9.6 


4 


07+° ii 


1.1+0.2 


92-393 


21: 


:39:44.08 


+57:42:15.9 


Low 


1451-1451-731-731 


M2.0 


2.0 


0.9 


9 


A 


1.0+0.2 


21395813 


21 


:39:58.13 


+57:28:33.5 


Low 


1 1 7-1 17-189-189 

_L_Li _L_LI A-fJJCf -i- O rj 


K2f 


1.7-'' 


8.0 


3 


A 


1.1+0.1 


93-720 


21: 


:40:09.99 


+58:00:03.6 


Low 


1451-1451-731-731 


m.Qf 


1.7/ 


1.8 


2 


A 


1.8+0.1 


24-1796 


21: 


:40:11.82 


+57:40:12.1 


Low 


1451-1451-731-731 


K7.0 


1.2 


9.5 


7 


"■"-0.3 


1.0+0.2 


21-2006 


21: 


:40;13.91 


+57:28:48.5 


Low 


44-44-63-63 


K5.0 


1.5 


12 


3 


^•^^-0.05 


1.7+0.3 


KUN-196 


21: 


:40:15.14 


+57:40:51.2 


Low 


122-122-144-144 


B9.0i 


2.2l 


4.0* 


D 




1.0+0.2 


21402192 


21: 


:40;21.92 


+57:30:05.4 


Low 


1451-1451-731-731 


K6.0 


1.4 


7.0 


9 


A 


1.5+0.1 


23-570 


21: 


:40:35.74 


+57:34:55.0 


Low 


1451-1451-731-731 


K6.0 


1.3 


6.4 


4 


^•'-0.3 


1.4+0.1 


23-162 


21: 


:40:44.52 


+57:31:31.9 


Low 


59-59-63-63 


K7.0 


1.7 


6.6 


1 


A 


1.9+0.1 


23-798 


21: 


:41;28.64 


+57:36:43.2 


Low 


975-853-377-377 


K5.5 


2.2 


15.0 


5 


0.79^ 


2.1+0.1 


NGC 7160 
























DG-481 


21: 


:52;21.10 


+62:45:03.7 


Low 


122-122-975-975 


A7.0 


1.5 


10.0* 


D 




1.6+0.1 


01-1152 


21: 


:53;19.82 


+62:34:00.4 


Low 


2418-2418-1950-1950 


Ml. 5 


0.9 


5.6 


10 


0* 


2.1+0.1 


01-580 


21: 


:53;37.07 


+62:28:47.2 


Low 


122-122-488-488 


K4.5 


1.6 


11.4 


6 


4+3 


2.7+0.1 



Comments 



Not detected 



SBl: 



SB2: 



CO 



Note. — Disked stars observed with IRS and their characteristics. Spectral types and extinctions from Contreras et al.(2002; marked with ^) and Sicilia-Aguilar et al. 
(2004, 2005); the data for 13-277 from Sicilia-Aguilar et al. (2008). Spectral types and extinction obtained from SED fitting (see text) are marked with /. The ages are 
derived from the V vs V-I diagram and the Siess et al. (2000) isochrones, except for the intermediate-mass stars (marked with *), for which we take it to be the cluster 
median. Early K-G stars and/or those affected by scattering in the V-I colors tend to have anomalously large ages, not taken into account in the analysis. The SED types 
refer to the model prototypes in Section 13.31 while intermediate- mass stars are marked as "D", given they are most likely debris disks or dust-depleted TD with large inner 
holes. Accretion rates and presence/ absence of accretion from Sicilia-Aguilar et al. (2005,2006b, 2010). Accretion rates marked with ' come from non-simultaneous data and 
may be more uncertain due to the intrinsic stellar variability. Objects marked with "A" have broad Ho and are thus accreting, although the accretion rate has not been 
measured. Objects marked with "PA" have strong Ho EW and are most likely accreting; conversely, objects marked with "0*" have small and narrow Ho emission and are 
probably not accreting. Uncertain values are marked with ":". The SED type number corresponds to the models in Table [TOl sec text. A few cases where the disks around 
intermediate- mass stars may correspond to debris disks are marked with "D". 
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Tablo 2. IRAM 1.3mm Observations 



Name 


Program(s) 


Nr. of scans 


Tau 


Flux (mjy) 


Mdiafc (M0) 


Comments 


Tr 37 














73-758 


132-07,145-08 


6 


0.2-0.4 


<1.8 


<0.009 




73-311 


132-07,145-08 


5 


0.2-0.3 


<1.8 


<0.009 




73-71 


132-07,145-08 


6 


0.2-0.3 


<1.5 


<0.008 




MVA-1312 


146-06 


3 


0.2-0.3 


<4.8 


<0.024 




KUN-314S 


146-06 


3 


0.2-0.7 


<6.3 


<0.032 


Nebulosity contamination 


11-2146 


146-06 


3 


0.2-0.7 


2.9±0.6 


0.015 




11-1209 


146-06 


3 


0.2-0.4 


<2.4 


<0.012 




11-2037 


146-06 


2 


0.2-0.3 


3.3±0.9 


0.017 




14-11 


146-06 


2 


0.2-0.5 


<2.4 


<0.012 


Nebulosity contamination 


21374275 


146-06 


1 


0.2 


<3.3 


<0.017 


Nebulosity contamination 


12-1091 


146-06,132-07,145-08 


7 


0.1-0.4 


<1.8 


<0.009 




MVA-426 


146-06 


2 


0.2-0.2 


2.2±0.7 


0.011 




82-272 


146-06,132-07,145-08 


3 


0.2-0.3 


<2.4 


<0.012 




13-669 


146-06,132-07,145-08 


3 


0.2-0.3 


<2.4 


<0.012 




13-1877 


146-06 


2 


0.3 


<2.7 


<0.014 




13-277 


110-06 


2 


0.4 


13.9±1.0 


0.070 


Sicilia-Aguilar et al. (2008) 


13-236 


146-06 


1 


0.3 


<3.9 


<0.020 




12-2113 


146-06,132-07,145-08 


4 


0.2-0.4 


<1.8 


<0.009 




13-157 


146-06,132-07,145-08 


4 


0.2-0.4 


<2.1 


<0.011 




54-1547 


146-06,132-07,145-08 


4 


0.2-0.4 


2.6±0.6 


0.013 




91-506 


132-07,145-08 


5 


0.2-0.5 


2.2±0.5 


0.011 




13-1250 


146-06 


2 


0.2-0.4 


2.4±0.8 


0.012 




21392541 


132-07,145-08 


3 


0.2-0.4 


<1.8 


<0.009 




24-515 


146-06,132-07,145-08 


6 


0.2-0.4 


<1.8 


<0.009 




21-33 


132-07,145-08 


5 


0.2-0.4 


<1.8 


<0.009 




92-393 


132-07,145-08 


5 


0.2-0.4 


<1.8 


<0.009 




21395813 


146-06,132-07,145-08 


3 


0.2-0.4 


<2.1 


<0.011 




93-720 


132-07,145-08 


4 


0.3-0.4 


<2.1 


<0.011 




KUN-196 


146-06 


2 


0.2-0.5 


<2.7 


<0.014 




21402192 


146-06,132-07,145-08 


3 


0.2-0.3 


<1.5 


<0.008 




23-798 


132-07,145-08 


4 


0.2-0.3 


<2.1 


<0.011 




NGC 7160 














DG-481 


146-06 


4 


0.4-0.6 


<2.4 


<0.012 




01-1152 


132-07,145-08 


7 


0.2-0.5 


<1.2 


<0.006 




01-580 


146-06 


4 


0.3-0.6 


<2.4 


<0.012 





Note. 



— Summary of IRAM 30m observations. To derive the total disk mass, the opacity of Beckwith et al. (1990; Ku 
and a gas to dust ratio of 100 have been applied. Upper limits correspond to 3(T for non-detections. 
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Table 3. Spitzer IRAC and MIPS fluxes 



Target IRAC 3.6/im (/ijy) IRAC 4.5/im (/ijy) IRAC 5.8/im (/ijy) IRAC 8.0/im (/ijy) MIPS 24/:im (/ijy) 



Tr 37 












73-758 


_ 


4264±63 


4308±206 


6811±91 


15638±101 


72-1427^ 


4663±50 


3721±54 


2932±143 


3013±93 


13831±484 


73-472-f 


9985±68 


8254±77 


5865±155 


3934±82 


24402±100 


73-311 


4373±48 


4451±66 


4205±156 


4388±88 


4418±101 


73-71 




8258±110 


6028±157 


5412±195 


9652±102 


72-875 


_ 








19814±1204 


21362507 


3060±58 


3013±61 


3021±224 


4872±26 


9196±103 


21364762 






_ 






14-141 


43079±170 


43062±137 


41906±335 


45029±172 




11-2146 


30102±130 


27881±122 


22736±266 


27283±152 


76509±112 


11-1209 


15601±68 


13757±82 


12274±176 


10310±110 


3950±104 


11-2322 


12399±82 


10551±79 


9468±223 


9556±39 


16711±106 


11-2037 


14656±90 


12443±61 


9938±245 


11975±82 


35664±106 


14-11 


7106±41 


5263±59 


3509±154 


3815±23 


26527±1881 


14-125 


10955±77 


10097±78 


9187±248 


9681±37 


28424±1388 


11-2131 


15234±91 


14013±65 


11430±243 


10879±82 


18036±482 


11-2031 


22307±112 


20279±106 


17741±266 


25481±153 


56282±933 


14-160 


11932±82 


10659±80 


9423±223 


11355±114 


26958±571 


14-183 


13177±83 


10027±78 


7779±212 


6151±103 


19010±104 


CCDM-l-5734 


249310±529 


196333±201 


106594±374 


95747±168 


107598±3171 


21374275 


57090±185 


44097±113 


50220±344 


64140±162 


129377±1549 


12-2519 


5857±64 


4796±63 


3945±208 


3818±91 


8514±103 


12-1091 


12677±84 


10811±80 


9228±240 


11157±104 


20119±103 


13-1238 


9730±62 


8990±61 


8593±207 


9645±91 


18905±101 


82-272 


32715±107 


28404±128 


23378±227 


26260±150 


75995±108 


13-669 


15991±93 


14291±67 


13162±245 


15883±89 


14891±103 


21380979^'^ 


3201±57 


2740±57 


2811±211 


3401±81 


11857±102 


54-1781 


3543±44 


2884±58 


2517±192 


1994±78 


3406±101 


13-1877 


17269±93 


15949±91 


14760±255 


16650±125 


35512±100 


13-277 


325705±757 


397223±456 


316238±775 


514053±552 


779325±2871 


13-236 


28167±131 


24697it85 


21195±250 


19573±107 


11950itl03 


12-2113 


19438±104 


19177±102 


20070±269 


24421±142 


34765±103 


13-157 


24698±117 


21381±83 


17888±244 


18270±130 


22978±104 


13-52 


3535±61 


2512±46 


2001±207 


1976±65 


4377±99 


91-155 


6372±66 


6618±58 


7301±241 


10165±86 


24939±100 


21384350 


2956±57 


2316±57 


2145±196 


1965±91 


2503±101 


54-1547 


8040±71 


6520±70 


5242±234 


5793±99 


16233±102 


91-506 


7893±68 


6689±69 


5267±249 


5561±89 


14074±97 


12-1617 


7328±52 


5780±54 


4475±168 


3983±68 


10447±100 


13-1250 


10045±77 


8116±72 


6939±217 


10332±108 


30610±105 


21392541^'*f 


12243±81 


9898±77 


7976±252 


9388±105 


14261±101 


21392570^'*^ 


4498±58 


3490±49 


2649±260 


2837±63 


9784±98 


21393104^'**^ 


21028±107 


16929±40 


12330±252 


14567±100 




24-515 


2889±43 


2121±58 


1763±153 


1705±87 


5134±99 


21-998 


8197±71 


7317±72 


6415±240 


5759±94 


10539±104 


21-33 


2407±57 


2266±48 


1956±213 


2096±74 


2711±97 


92-393^*^ 


4097±58 


3371it61 


2446±262 


2589±77 


6408±144 


21395813^'*^ 


38709±157 


33078±142 


27382±269 


23714±143 


16536±101 
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Table 3 — Continued 



Target IRAC 3.6/im (/ijy) IRAC 4.5/im (/ijy) IRAC 5.8/im (/ijy) IRAC 8.0|Um (/ijy) MIPS 24;im (juJy) 



93-720 












24-1796 


3995±58 


2945±53 


2409±221 


1850±70 




21-2006 


7542±60 


6835±54 


5943±234 


7128±74 




KUN-196 


5455±63 


3560±61 


2214±214 


1476±85 




21402192 


3801±57 


2877±58 


2671±218 


2482±82 


2648±100 


23-570 


7886±69 


6878±70 


6178±248 


6064±94 




23-162 


7328±67 


6479±51 


6260±260 


6430±62 


9797±141 


23-798 












NGC7160 












01-1152 


1113±51 


1083±52 


1069±244 


1424±67 


928±142 


01-580 


8255±13 


7319±62 


5530±146 


4862±67 


8125±307 


DG-481 


9813±65 


7671±71 


6574±207 


5376±82 


2395±187 



Note. — IRAC and MIPS fluxes for the targets observed with IRS. Most of the objects had aheady been presented in SA06, 

and their newly reduced values are in good agreement with those. New IRAC and MIPS detections (due to the inclusion of 
archive data from the region, better data from the new pipelines, and/or new members found in Sicilia-Aguilar et al. 2006b) 
are marked with ^ and , respectively. IRAC and MIPS photometry for 21364762 is contaminated by the strong ba<;kground, 
and 93-720 is located out of the IRAC and MIPS fields. 
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Table 4. Optical data 



Target U V R I References 



Cousins System photometry 
Tr 37 



73-758 




17.062 


15.845 


14.761 


1 


72-1427 





18.080 


16.822 


15.580 


1 


73-472 





16.851 


15.755 


14.805 


1 


73-311 


_ 


18.378 


17.126 


15.887 


1 


73-71 




16.989 


15.813 


14.694 


1 


72-875 




18.367 


17.382 


16.178 


1 


21362507 


_ 


20.426 






1 


14-141 


19.349 


15.810 


14.680 


13.519 


1 


11-2146 


18.157 


16.917 


15.688 


14.412 


1 


11-1209 


17.695 


15.414 


14.522 


13.596 


1 


11-2322 


18.690 


16.687 


15.591 


14.383 


1 


11-2037 


18.115 


16.033 


15.085 


14.127 


1 


14-11 




17.528 


16.145 


14.774 


1 


14-125 


20.027 


16.738 


15.733 


14.701 


1 


11-2131 


19.815 


17.722 


16.433 


15.248 


1 


11-2031 


18.084 


15.483 


14.582 


13.697 


1 


14-160 


20.234 


16.909 


15.843 


14.832 


1 


14-183 


20.838 


16.613 


15.371 


14.217 


1 


CCDM-l-5734 




8.240 







3 


12-2519 




17.274 


16.231 


15.210 


1 


12-1091 


18.613 


16.587 


15.652 


14.730 


1 


13-1238 


19.149 


18.355 


16.827 


15.370 


1 


82-272 


_ 


17.257 


16.170 


14.997 


1 


13-669 


18.285 


15.836 


14.857 


13.922 


1 


21380979 




18.793 






1 


54-1781 


_ 


18.180 


16.999 


15.654 


1 


13-1877 


19.190 


16.829 


15.591 


14.445 


1 


13-277 


16.987 


13.915 


12.952 


12.141 


1 


13-236 


17.559 


15.641 


14.706 


13.856 


1 


12-2113 


18.786 


17.143 


15.920 


14.669 


1 


13-157 


17.510 


16.234 


15.255 


14.327 


1 


13-52 




17.265 


16.182 


15.148 


1 


91-155 




18.306 


16.943 


15.515 


1 


21384350 




20.025 






1 


54-1547 


18.655 


16.653 


15.698 


14.813 


1 


91-506 




16.842 


15.766 


14.756 


1 


12-1617 


19.721 


17.786 


16.508 


15.189 


1 


13-1250 


18.676 


16.049 


15.090 


14.235 


1 


24-515 




17.896 


16.700 


15.576 


1 


21-998 


19.178 


17.478 


16.395 


15.301 


1 


21-33 




18.677 


17.397 


16.141 


1 


92-393 




18.272 


16.826 


15.354 


1 


93-720 




16.453 


15.368 


14.362 


1 


24-1796 


19.157 


17.395 


16.367 


15.312 


1 


21-2006 


19.661 


16.994 


16.005 


15.020 


1 


KUN-196 




13.850 






3 


23-570 




16.847 


15.715 


14.873 


1 
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Table 4 — Continued 



Target 


U 


V 


R 


I 


References 


23-162 


— 


17.974 


16.705 


15.689 


1 


23-798 


19.717 


18.755 


17.423 


16.360 


1 


NGC 7160 












01-1152 


— 


18.574 


17.370 


16.186 


1 


01-580 


— 


16.937 


15.958 


15.054 


1 


DG-481 


— 


13.220 


113.220 


113.220 


2 


Johnsons System Photometry 












Tr 37 












21364762 


— 


18.472 


16.809 


15.116 


4 


21374275 


— 


17.485 


16.355 


15.128 


4 


21392541 


— 


16.577 


— 


13.733 


4 


21393104 


— 


15.690 


— 


13.326 


4 


21395813 


— 


14.545 


— 


12.664 


4 


21402192 


— 


16.964 


15.415 


14.485 


4 


72-1427 


20.176 


18.272 


16.618 


15.161 


4 


72-875 


20.814 


18.245 


16.590 


15.006 


4 


21362507 


21.833 


20.426 


18.578 


— 


4 


14^141 


19.219 


16.036 


14.508 


13.094 


4 


11-2146 


18.684 


16.662 


15.211 


13.999 


4 


11-1209 


17.755 


15.340 


14.040 


13.094 


4 


11-2322 


18.858 


— 


15.357 


13.953 


4 


11-2037 


18.144 


— 


— 


13.687 


4 


14-11 


20.899 


17.653 


16.170 


14.269 


4 


14^125 


19.169 


16.616 


15.224 


14.041 


4 


11-2131 


20.123 


— 


15.853 


14.605 


4 


11-2031 


17.601 


— 


14.240 


13.345 


4 


14-183 


19.838 


17.112 


15.648 


14.431 


4 


12-2519 


20.539 


— 


15.958 


14.763 


4 


12-1091 


18.983 


16.274 


14.984 


13.984 


4 


13-1238 


19.259 


17.902 


16.337 


14.601 


4 


82-272 


19.414 


16.616 


15.266 


14.019 


4 


13-669 


18.231 


15.834 


14.630 


13.507 


4 


54-1781 


20.299 


18.564 


16.922 


15.120 


4 


13-1877 


18.883 


17.189 


15.620 


14.222 


4 


13-277 


15.016 


13.146 


12.090 


11.117 


4 


13-236 


17.665 


15.546 


14.431 


13.412 


4 


13-157 


17.861 


16.768 


15.320 


14.012 


4 


21384350 


23.366 


19.918 






4 


54-1547 


18.539 


16.571 


15.299 


14.437 


4 


91-506 


19.777 


17.141 


15.732 


14.546 


4 


12-1617 


20.831 


18.039 






4 


13-1250 


18.403 


15.956 




13.668 


4 


24^515 


19.908 


17.941 


16.396 


15.123 


4 


21-998 


19.798 


17.668 




14.961 


4 


21-33 


20.814 


18.723 


17.138 


15.722 


4 


24-1796 


18.389 


17.341 


15.897 


14.825 


4 


21-2006 


19.564 


16.933 


15.517 


14.589 


4 


23-570 


18.391 


16.718 


15.276 


14.308 


4 
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Table 4 — Continued 



Target 


U 


V 


R 


I References 


NGC 7160 










01-580 


18.380 


17.392 


15.755 


14.737 4 



Note. — Suiiiiiiary of optical magnitudes for the analyzed sources. 
References: (1) Sicilia-Aguilar et al. (2005); (2) De Graeve (1983); 
(3) Marschall & van Altena(1987); (4) Sicilia-Aguilar et al. (2010). 
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Table 5. 2MASS data 



Target 



2MASS ID 



H 



K 



Tr 37 

73-758 

72- 1427 

73- 472 
73-311 
73-71 
72-875 
21362507 
21364762 
14-141 
11-2146 
11-1209 
11-2322 
11-2037 
14^11 
14-125 
11-2131 

11- 2031 
14-160 
14-183 

CCDM-l-5734 
21374275 

12- 2519 

12- 1091 

13- 1238 
82-272 
13-669 
21380979 
54-1781 
13-1877 
13-277 
13-236 

12- 2113 

13- 157 
13-52 
91-155 
21384350 
54-1547 

91- 506 

12- 1617 

13- 1250 
21392541 
21392570 
21393104 
24^515 
21-998 
21-33 

92- 393 
21395813 



21350835-1-5736028 
21351627-1-5728222 
21351861+5734092 
21352451-1-5733011 
21353021-1-5731164 
21354975-1-5724041 
21362507-1-5727502 
21364762-1-5729540 
21364941+5731220 
21365767+5727331 
21365850+5723257 
21370191+5728222 
21370703+5727007 
21371031+5730189 
21371054+5731124 
21371215+5727262 
21371591+5726591 
21372732+5731295 
21373849+5731408 
21374093+5733372 
21374275+5733250 
21375107+5727502 
21375762+5722476 
21375926+5736162 
21380350+5741349 
21380928+5733262 
21380979+5729428 
21381612+5719357 
21381703+5739265 
21381731+5731220 
21382742+5731081 
21382743+5727207 
21382804+5730464 
21383255+5730161 
21383470+5741274 
21384350+5727270 
21384446+5718091 
21385807+5743343 
21390468+5725128 
21391213+5736164 
21392541+5733202 
21392570+5729455 
21393104+5747140 
21393407+5733316 
21393480+5723277 
21393561+5718220 
21394408+5742159 
21395813+5728335 



13.330+0.028 
14.005+0.031 
13.219+0.024 
14.014+0.034 
12.979+0.034 
14.092+0.031 
14.952+0.044 

13.568: 
11.918+0.022 
12.442+0.031 
12.200+0.024 
12.854+0.029 
12.649+0.026 
13.057+0.029 
13.090+0.026 
13.393+0.027 
12.499+0.025 
13.218+0.037 
13.303+0.029 
7.697+0.019 
12.561+0.023 
13.647+0.030 
13.182+0.024 
13.321+0.027 
12.502+0.026 
12.391+0.029 
14.161+0.040 
13.991+0.032 
12.935+0.039 
10.279+0.026 
12.355+0.036 
12.856+0.030 
12.912+0.030 
13.738+0.047 
13.741+0.035 
14.246+0.034 
13.421+0.027 
13.386+0.031 
13.542+0.026 
12.920+0.024 
12.590+0.025 
13.942+0.032 
12.410+0.027 
14.057+0.028 
13.850+0.029 
14.501+0.039 
13.619+0.026 
11.720+0.024 



12.533+0.032 
13.067+0.029 
12.299+0.031 
13.115+0.037 
12.105+0.037 
13.006+0.027 
14.031+0.048 
12.342+0.042 
10.914+0.029 
11.327+0.032 
11.433+0.033 
11.969+0.033 
11.752+0.027 
12.216+0.034 
12.170+0.029 
12.214+0.031 
11.534+0.028 
12.349+0.041 
12.229+0.031 
7.647+0.031 
11.308+0.026 
12.816+0.036 
12.215+0.029 
12.400+0.032 
11.483+0.029 
11.566+0.029 
13.356+0.044 
13.092+0.032 
11.954+0.043 
9.329+0.029 
11.363+0.040 
11.986+0.035 
11.979+0.036 
12.842+0.040 
12.875+0.048 
13.447+0.037 
12.647+0.035 
12.521+0.036 
12.627+0.028 
12.119+0.029 
11.686+0.029 
12.996+0.032 
11.491+0.028 
13.125+0.033 
12.972+0.033 
13.581+0.040 
12.765+0.028 
10.968+0.031 



12.163+0.028 
12.673+0.028 
11.829+0.022 
12.674+0.030 
11.715+0.029 
12.500+0.026 
13.518+0.040 
11.655+0.033 
10.355+0.021 
10.640+0.025 
11.122+0.023 
11.560+0.025 
11.258+0.020 
11.956+0.026 
11.732+0.021 
11.523+0.025 
10.856+0.020 
11.954+0.033 
11.668+0.045 
7.611+0.021 
10.388+0.021 
12.460+0.022 
11.673+0.020 
11.889+0.018 
10.847+0.019 
11.195+0.020 
12.974+0.042 
12.792+0.030 
11.263+0.024 
8.593+0.020 
10.774+0.024 
11.506+0.024 
11.249+0.052 
12.607+0.031 
12.489+0.036 
13.065+0.034 
12.180+0.026 
12.040+0.022 
12.129+0.024 
11.799+0.023 
11.312+0.018 
12.700+0.028 
11.023+0.023 
12.704+0.028 
12.476+0.028 
13.220+0.044 
12.505+0.023 
10.422+0.022 
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Table 5 — Continued 



Target 


2MASS ID 


J 


H 


K 


93-720 


21400999+5800036 


12.709+0.024 


11.580+0.030 


10.940+0.022 


24-1796 


21401182+5740121 


13.876+0.035 


13.076+0.039 


12.734+0.037 


21-2006 


21401390+5728481 


13.704+0.029 


12.920+0.038 


12.488+0.030 


KUN-196 


21401508+5740513 


12.145+0.035 


11.922+0.046 


11.881+0.051 


21402192 


21402192+5730054 


13.608+0.039 


12.754+0.038 


12.538+0.032 


23-570 


21403574+5734550 


13.477+0.030 


12.556+0.030 


12.194+0.028 


23-162 


21404450+5731314 


14.044+0.030 


13.144+0.037 


12.588+0.029 


23-798 


21412864+5736432 


14.410+0.022 


13.533+0.035 


12.820+0.030 


NGC 7160 










01-1152 


21531982+6234004 


14.930+0.044 


14.078+0.050 


13.764+0.043 


01-580 


21533707+6228469 


13.707+0.029 


12.887+0.031 


12.301+0.025 


DG-481 


21522113+6245034 


11.918+0.027 


11.713+0.031 


11.520+0.021 



Note. — Summary of 2MASS magnitudes for the analyzed sources. Most of the data 
are already published in Sicilia-Aguilar et al. (2004, 2005). 
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Table 6. Silicate Composition in the 7-14/xm region 



Name 


O.l/im A.Ol. 


0.1/im A.Py. 


0.1/im Fors. 


0.1/im Enst. 


0.1/im Sil. 




1.5;im A.Ol. 


1.5/im A.Py. 


1.5//m Fors. 


1.5/im Enst. 


1.5/(m Sil. 




6.0;im A.Ol. 


6.0/im A.Py. 


6.0//m Fors. 


— 


6.0/im Sil. 


TV. QT 












73-758 


"■'-0.7 




2 0+°^'' 
^■"-0.3 


1+°-* 

"■-'-0.1 






12+5 
-6 


17+11 

'-8 


7+^? 
"■ ' -0.6 


1 6+;^? 


2+°^I 
" -0.2 




7+18 
'-7 


55+1° 






4+4 
^-3 


73-472 


30+21 

'^"— 20 


9+12 


2 5+;^3 

^-1.2 


5+i-2 


1+°-? 

" -0.1 




91+55 


8+9 

°-7 


6+16 

"■"-0.6 


9 7+3.4 
''•'-2.3 


4+1-2 
"■^-0.4 




9K+38 








^■1-0.1 


73-311 


2+3 


3+5 




1 5+0-« 

^•"-0.8 


5+°6 
■ -0.3 




2+2-9 


19+1° 


2+11 
"■''-0.2 


1+2-2 
"•^-0.1 


1+13 
"•^-0.1 




48+16 


8+25 
°-8 


3+8 




15+« 


73-71 




11+2 


6+°-2 
"■"-0.2 




1 3+03 




2+?-2 

" -0.2 


8+1! 

"-4 


3+2-6 

—0.3 


3 9+0-? 

—0.9 






71 +4 


r\ 0+5.3 
"■'*-0.3 


q+2.6 
0.0_2.3 




n 9+0-'' 

"•^-0.2 


21362507 


17+i 

' -5 


4+^ 
-3 










3+^ 

^-2 


10+ii 
"-7 


5+2 

^-2 


2+1 
-1 








5»_16 


9+14 
^-2 






11-2146 


12+3 


1+4 
-'■-I 


1 4+0-6 
-^■^-0.4 




"■' -0.4 




16+23 


99+10 
^^-10 




3+2 

•J-2 






11+16 


32+1* 


<9 




<4 


1 1 9n'?7 

1 ±-ZUO ( 


<.D 






9+1 
^-1 


0+1 




1+9 
1-1 


42+19 


1+2 
1-1 








09+21 


16+2° 






1+3 
^-1 


14-125 




10+6 
1U_3 


9+°6 
"■^-0.4 


8+°-'' 
"•°-0.5 


1 5+°-9 
■l-^-0.5 






2+'' 
^-2 


2+11 
"-^-0.2 


4+2 
^-1 






15+" 


64+1° 






2+'* 
''-I 


11-2031 


16+3 
1D_3 


<3 


2+1 
^-1 








<12 


38+12 


7+2 
' -2 




2+1 
''-I 




12+1* 


16+15 

-'■"-13 






4+6 
^-4 


14-160 


32+20 




8+°^4 

"■"-0.3 








5+'! 


6+16 

"-6 


3+1 








49+18 






4+11 
^-4 


14-183 


68+" 




4+1 
^-1 










4+9 


4+3 

— 2, 








<35 










21374275 






l-=-0.2 




0.2t°;l 






14+1 


2 0+°^5 
^■"-0.6 


1 r+0.5 
1-^-0.4 


6+°3 
"•"-0.3 




0.813-6 


79+1 








12-2519 


0.2t6;3 


9+8 
^-2 


0+4 
•5-2 


1+1-2 
"•-'■-0.1 


2+'' 
^-1 






1+14 
■'-I 


4+10 
^-3 


0+10 
*»-5 


4+11 
^-3 




QQ + 25 


^^-20 


9+9 
^-2 




15+22 

i^-io 


12-1091 




9+10 


411 








]^+15 


24+13 


It? 


613 


9+2 






1+21 
■'-I 






lit 


13-1238 
















90+37 


It? 


7+3 
'-2 


0+5 

•1-2 






50+2'' 









Table 6 — Continued 



Name 0.1/^m A.Ol. 0.1/im A.Py. 0.1/im Fors. 0.1/im Enst. 0.1/im Sil. 

1.5/im A.Ol. 1.5/im A.Py. 1.5/im Fors. 1.5/im Enst. 1.5/im Sil. 
e.O^im A.Ol. 6.0/im A.Py. 6.0/im Fors. — 6.0//m Sil. 



82-272 15 



g+18 

91-506 62+^ 



4 



+8 ,+18 Q+1 9+2 9+1 

-7 "^-3 °-l ^-2 ^-1 

<11 60tfi — 

<3 — — 

13-669 — — 2.4+Q2 0.7lo4 l-2loi 

44^2 1+3 

32380979 2lf — — — 

54111 - 1511 
90+I6 

13-1877 0.ll°? 0.81^-2 2.5t°| - 

— ' 49+^' 1+°''' 4+^' 0&+°-^ 

^=■-9 ^-^-0.1 ^-1 ^-"-O-S 

9n~t~^^ 9n~t~^8 9+2 

^Ll-ll ^"-15 ^-2 

13-277 23^1 — l-8lo2 O-ltoi ~ 

— 43+* 1 1+°-'^ 2 l+°-'' 1 0+''-* 

^•J_3 -L-J^_o.7 ^-l-O.T ^-"-0.4 

_ 27+4 _ _ 4+^'^ 

13-236 — — 0.7l°l O.ltoi O.Stol 

— 3 5+°-^ 1 4+°-* 1 fi+° '* 
•5-0_o.8 -'■■^-0.4 '■■"-OA '^•-'^-0.1 

57tl 34l^ _ _ _ 

12- 2113 34^^ <4 — — — 

2li-' 36t\l - ll? - 

Q+26 1fi+19 1+5 

»_9 ^°-14 ^-1 

13- 157 22^^ 2tl 4tl — — 

q+20 qn+13 9+2 i;+4 

0_3 •^"-13 ^-2 ^-3 

c+38 9+27 9^+6 

13-52 — <3 3^1 3ll — 

16l§ - 5il 

66l?o 711* _ _ _ 

91-155 40II2 <3 3ll — — 

7lf 2911^ - 131^ _ 

— itT <5 - btl 
54-1547 66lg itl 4tJ 2^1 — 

— letl - it? 3ll 

<5 — — <4 

2ll^ _ _ _ 

4+7 0+15 p.+2 

^-4 '^-3 "-2 

12- 1617 4^3 6^4 3^1 2t\ — 

<4 49i;° - 5tl - 

13+\l _ _ _ iv+4 

13- 1250 2ll^ — 3lJ — — 

fi+16 Q1+14 9+1 p;+2 

°-6 '^^-9 ^-1 "^-2 

9+12 'jn+ll 1+8 

^-2 '^"-20 ^-1 

21392541 24^2 



38 



+8 q+1 



20l^ lOl^ - - 2tl 
+ 10 9+19 ,+1 



21392570 2+i" 2 



2 ^-2 "^-1 



Table 6 — Continued 



Name 0.1/im A.Ol. 0.1/im A.Py. O.l^tm Fors. 0.1//m Enst. O.l^tm Sil. 

1.5/im A.Ol. 1.5jUm A.Py. 1.5/im Fors. 1.5/im Enst. 1.5/im Sil. 
6.0/im A.Ol. e.O^im A.Py. 6.0;im Fors. — 6.0;im Sil. 



21393104 



21-' 



21395813 



93-720 



21402192 



23-570 



23-162 



23-798 



NGC 7160 

01-1152 



01-580 



DG-481 



11-1-6 
°^-31 

4711 



o.i; 



o.i; 



-1.3 
-0.1 



-1.1 
-0.1 



10-H6 
^°-14 



+21 



21 



7tl 



+32 
-1 
I -1-11 

'--14 



10 
19 



+41 
-10 
+2 
2 



+8 



52 



<3 
4lf 
43l^ 

q-l-12 

Sir 

8tl' 



16 

6 
+ 24 
-1 
1.1 
0.1 

4411 
<4 



12 



1 

0.1 



91 



+3 
5 



17+27 
'-'-15 

0+22 

•^-3 



35l^ 
2011 



91 



65 



It? 

+3 
-6 
10 
4 
+8 

15 
+43 
-3 



39 
3 



+10 

16 
+29 
3 



13 



+s 



□+9 



hll 
-16 
5+23 
-13 
,1+10 

hlO 
-9 
hl2 



4! 
49l 
13^ 



16^ 
6j 



31^ 

Q 9+0.2 
•^•^-0.2 



7+0.3 
^- ' -0.3 



+3 

3 



+6 
4 



5 

3li 

1+3 



11 

0.5 
1 



1 

+8 
-7 
0.1 
0.1 
+0.2 
-0.2 



+ 1 
-1 
+3 
-2 



1.3 



+0.2 
0.2 

211 



56 



+7 



2II 
4l^ 
<7 

3li 
2II 
ll? 



2+°-^ 
'•'■■^-0.2 

n f,+0.8 



6+°" 



121^ 



4+1 



41? 



2II 
It? 



2II 



.+2 



2II 



311 



6ir 



1 q+0.3 



'^••^-0.3 
4+0-3 



6l^ 
21^ 



6 3+° ''' 

U.D_Q 2 

o.5l°o;! 
It? 

711 
211 

It? 



It? 



^-2 



16 



+6 
6 



7+3 
1-4 



7+4 
'-4 



Note. — Silicate fit for the objects in Ccp OB2 using the Two Layer Temperature 
Distribution metfiod (J09) in the short wavclengtfi range (7-14/im). Only those objects 
with enough S/N to significantly constrain the dust mineralogy are listed here. The values 
are percentages in mass fraction. If no value is given for a certain species/size, it means 
that it is present at 2% level or less. 



Table 7. Silicate Composition in the 17-35//m region 



Name O.l/um A.Ol. 0.1/xm A.Py. 0.1/xm Fors. 0.1/xm Enst. 0.1/xm Sil. 

1.5/xm A.Ol. 1.5/im A.Py. 1.5;um Fors. 1.5/xm Enst. 1.5/xm Sil. 
6.0/xm A.Ol. 6.0fim A.Py. 6.0/im Fors. — G.Ofim Sil. 



Tr 37 

11-2146 6+f 10t^° Si? — — 

49l?8 <7 2li — 

— 25tll <7 — ^-l 



CCDM+5734 — 19^}^ 5tJ 

6lt? lll{? 

9+6 
^ -2 



13-277 



82-272 — 54:^11 8^1 2^1 



Utt 18111 - 

O+60 9+42 9+24 

0_2 



621:1 181:1 



18tl - - - Itl 

21393104 471} — 31} — — 

- 441^ — 31} - 

- <4 — — — 

93-720 42l}^ 20l[i ^1} — — 

261^2 3+f _ _ _ 

- 311 - - li? 



Note. — Silicate fit for the objects in CepOB2 using the Two Layer Temperature 
Distribution method (J09) in the long wavelength range (17-35;um). Only objects with 
good S/N in this region are presented. The values are percentages in mass fraction. If 
no value is given for a certain species/size, it means that it is present at 2% level or less. 



Table 8. Grain sizes and crystalline fraction in the 7-14/im region 



Name 



Size (am.) / fj.m Size (cryst.) / fim Cryst.Fract. Forst/Enst 



Comments 



Tr 37 








73-758 




n fi+0.1 
0-8-0.4 


A+2 
4-1 


TQ A '70 

/ 0-4 / Z 


1 7-1-1.9 


n 7+0.2 
0-7-0.4 


c+S 


73-311 


A -7+0.3 


O-elas 


Q+5 

3Ii 


73-71 


A n+0-2 


q 9 + 0.6 

■^-^-1.3 


0+3 


zldozouY 


/I Q-1-0.3 


1 c + l.O 
^-3-0. 2 


•7+6 


11-2146 


q R + 0.3 


1 1+0.2 
1-1-0.4 


4-2 


11-zUo / 


3-2 


1 n+0.1 
1-0-0.1 


19+6 


14-125 


c c:-f0.2 

■^•■^-o.e 


1 1+0.2 
1-1-0.2 


5-1 


11-2031 


g r,-|-0.5 

3-0-1.2 


1 n+0.1 

1-2-0.1 


10+2 


14-160 


A 1-1-0.3 


1 1+0.1 
1-1-0.2 


Q+2 

3Ii 


14-183 


n /I+2.1 


n s+0.2 
0-8-0.4 


8-3 


ilv5/4z 10 


C Q + O.l 

''-'^-0.1 


1 n+Ol 
1-0-0.1 


c ^+0.7 

5-4-0.7 


12-2519 


6.0-2.7* 


1 q + l.O 

1-3-0.1 


0+20 

9-3 


12-1091 


3-9-1.0 


1 n+0.1 
1-0-0.2 


11+5 
11-3 


1 Q 1 OQQ 

lo-lzoo 


5 3+°-^ 


n o+O.l 
0-9-0.2 


1/1 + 10 
14-4 


or) 0*70 


1 n + Oil 
1-2-0.2 


n K+01 


1K+3 

lb 3 


io-Doy 


q Q-(-0.1 


1 1+0-7 
1-1-0.1 


I- 1+2 

II- 2 


21380979 


n Q + 0.7 


1 n+0.1 
1-0-0.1 


91 +7 

21-5 


13-1877 


q K+0.5 
''-^-O.S 


n a+0.1 
0-8-0.1 


Q + 1 
8I2 




9 c+0.2 
^-^-0.2 


1 n+0-1 

1-0-0.1 


5-1 


13-236 


0-1-0.2* 


1 9+0.1 
1-2-0.1 


q n+0-5 

3-9Io.6 


12-2113 


O q-l-0.9 

2-3-1.1 


n 9+0.8 
0-2-0.1 


1+2 
1-1 


13-157 


n n+10 
0-9-0.2 


1 n+0.2 
1-0-0.3 


in+5 
IOI4 


13-52 


c q+0.2 


7+°-l 


10+3 


91-155 




1 9+0.2 
^-^^-0.2 


16l^ 


54-1547 


o-5lS:I 


4+°-^ 




91-506 


6+°-^ 
^-°-o.i 


0.51°- 


5 9+1-0 


12-1617 


1 5+1-2 


0.810:1 


101^3 


13-1250 


q 1+0.3 
''-^-l.O 


1 l+°-i 




21392541 


2 g+0.3 


1 l+°-i 
^■^-0.1 


511 


21392570 


c: 9+0.3 


i-3lS;i 


181^ 


21393104 


o.8i°:i 


0.8lO;i 


6-4ir9 


21-998 


6.0* 


0+1 


stl 


21395813 


q c: + 0.9 


q q+0.6 

'-'■'-'-1.2 


24111 


93-720 


A 0+0.1 


1 2+°-l 


9 q + 0.2 
^-'-'-0.2 


21402192 


6.0-0.1* 


1 9+0.2 
^-^-0.1 


4.9irs 


23-570 


1 5+°-* 


6+°-2 
"■"-0.3 


16ll 


23-162 


9+1-5 

U.»_Q 2 


Oj+0.1* 


5±^2 


23-798 


^-^-0.4 


1 l+°-l 

^■^-0.2 


5il 



+2 

-2 
+4 
-2 
+4.3 



0-2-0.2 

1 0+"" 

q+5 

'^-3 

7+1-3 
' -0.7 

5+0-4 
^■^-0.1 

0+12 

Q 0+0.3 
q4+320 

311 

1+4 

7+0.8 
-0.4 

6+1-° 
^■°-0.4 

^-'-0.2 

5+0-3 
U.O_g 2 

5+0-2 
U.O_Q 2 

j+0.8 

-0.8 

1+0.5 



1.0! 
0.7^ 



1-3! 

1-2-0.5 

41^^ 

It? 
4+°-3 

,+0.3 
-0.2 

2+1 
^-1 
=+0.4 
'-0.3 
;+0.3 
-0.2 
=+0.7 
'-0.5 
,+0.3 
-0.2 
1+0.4 
-0.2 
-1+0.3 
-0.3 
1+5 
-3 

111 



0.3"' 



O.i 
0.5! 
O.i 
0.6! 
0.4j 
1.0^ 
4-* 



0.2 
1 
2 

1.5 



+0.1 
0.1 

1 

1 

+2 
2 

+0.9 
0.8 



TDn 

TDa, Low S/N 



Variable M, Continuum uncertain 

Variable M 
Low S/N 

Continuum uncertain 

Variable M 
Variable M 

TDa; 

Variable M 

Variable M, Continuum uncertain 



TDn 
Low S/N 



Variable M 

TDa, Continuum uncertain 
TDa; 

Continuum uncertain 
Continuum uncertain 

TDa, Continuum uncertain 



Continuum uncertain 



NGC 7160 

DG-481 
01-1152 
01-580 



4.7 
1.6 
0.9 



0.2 
0.2 

+ 1.7 
0.2 
+0.9 
0.3 



1.1 
1.2 

0.7 



+0.1 
-0.2 
+0.2 
0.1 
+0.8 
0.3 



7.3 
11 
7" 



+1.9 
-1.7 

+4 
3 
+4 
2 



1.3 
1.2 
2.3 



+0.8 
0.8 
1.0 
0.8 
3.0 
1.7 



TD/D, Continuum uncertain 
TDn; 



Note. — Average grain sizes for amorphous and crystalline silicates and crystalline fraction (mass fraction of crystalline 
silicates) from the TLTD fit in the 7-14/xm region. Only those objects with enough S/N to significantly constrain the dust 
mineralogy are listed here. The colon (;) indicates uncertain values in spectra with poor S/N. Grain sizes of 0.1 or 6/im 
indicate the presence of very small or large grains, but given that our model does not include smaller nor larger grains, the 
uncertainties cannot be well determined, we thus mark these objects with The comments include whether the star 

presents a variable accretion rate ("Variable M), whether it is an accreting or non-accreting transition disk (TDa, TDn, 
respectively), and whether the fitted continuum level may be uncertain as described in Section 13.31 
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Table 9. Grain sizes and crystalline fraction in the 17-35/um region 



Name Size (am.) / //m Sizo (cryst.) / fan Crysl. ( Forst/Enst) Commonts 



Tr 37 



11-2146 


9 c+0.6 


0.81°:? 


5l? 


1 q + 2.8 

-^■"-•-o.e 


CCDM+5734 


1 O+0.2 


0.1* 


c 1 +0.4 
^•-■--0.5 




82-272 


U.D_Q 2 


Q -|^+0.4* 


8 q+^-^ 

°-^-0.9 


3+? 


13-277 


9 q+0.3 


4+°-2 


1.01°:} 


2+? 


21393104 


^-3 




51? 


21^ 


93-720 


^•'-0.3 


4+°-2 


4.81?:? 


1 q+0.6 
-^•"^-0.6 



Note. — Average grain sizes for amorphous and crystalline silicates and crystalline fraction 
(mass fraction of crystalline silicates) from the TLTD fit in the 17-35;um range. Only objects 
with enough S/N in this region are presented. The intermediate-mass star CCDM-l-5734 is 
marked with TD/D since it is probably a transitional disk with a very large inner hole or a 
debris disk. The colon (:) indicates uncertain values in spectra with poor S/N. Grain sizes of 
0.1 jum indicate the presence of very small grains, but given that our model does not include 
smaller grains, the uncertainties cannot be well determined, we thus mark these objects with 



Table 10. Summary of stellar and SED properties of the disk prototypes 



Model 


Prototype 


Te// (K) 


R. (Ro) 


M* (Mq) 


Mdisk (Mq) 


Comments 


Normal (1) 


11-2037 


4770. 


2.3 


1.3 


0.015 


Normal disk with silicates 


Massive (2) 


13-277 


5520. 


5.9 


2.6 


0.07 


Very massive, strongly accreting disk 


Normal/Depleted (3) 


13-236 


4900. 


2.3 


1.6 


<0.020 


Settled disk with silicates 


Settled/Depleted (4) 


11-1209 


4205. 


2.9 


1.0 


<0.012 


Settled disk with no/marginal silicates 


Kink (5) 


13-1250 


4470. 


2.2 


1.3 


0.012 


Disk with strong slope change or kink 


Old CTTS disk(6) 


01-580 


4470. 


1.5 


1.2 


<0.012 


Depleted disk with strong silicate feature 


Classical TD (7) 


24-1796 


4060. 


1.4 


0.9 




Accreting TD without silicate, strong turn-up 


Classical TD, sil (8) 


14-11 


3650. 


2.1 


0.4 


<0.001 


Non-accreting TD with silicate, strong turn-up 


Settled TD, no sil (9) 


92-393 


3580. 


2.1 


0.6 


<0.009 


Accreting TD without silicate, settled 


Settled TD, sil (10) 


13-52 


4060. 


1.5 


0.8 




Non-accreting TD with silicate, settled 



Note. — Modeled disk prototypes. Temperatures use the correspondence between spectral types and Tgff from Kenyon & Hartmann (1995). 
Spectral types, stellar radii and masses from Sicilia-Aguilar et al. (2006b, 2010). 



Table 11. SED slope at different wavelengths 



Name a{H-5.8) a(H-8) a(H-12) a(H-17) a(3.6-5.8) a(H-24) a(5.8-8) a(8-17) a(24-33) 



73-758 


-1.77 


-1, 


.47 


-1, 


.16 


-1, 


.01 






-0, 


.94 


-0.26 


-0.05 


-1.78 


72-1427 


-1.81 


-1, 


.79 


-1, 


.38 






-1, 


.93 






-1.71 






73-472 


-1.85 


-1 


.98 


-1 


.78 


-1 


.81 


-2, 


.08 


-0, 


,95 


-2.50 


-1.45 


0.36 


73-311 


-1.50 


-1 


.42 


-1 


.21 


-1 


.60 


-1, 


.05 


-1, 


,46 


-1.10 


-2.00 


2.79 


73-71 


-1.79 


-1 


.84 


-1 


.62 


-1, 


.47 






-1, 


,32 


-2.04 


-0.68 


-1.75 


72-875 


-1.95 


-1, 


.35 


-1, 


.25 


-0, 


.80 






-0, 


,78 


0.99 


0.37 


0.42 


21362507 


-1.08 


-0 


.88 


-0 


.71 






-1, 


.06 






-0.08 




_ 


21364762 


-1.33 


-1, 


.30 


-1, 


.76 










-0, 


.57 


-1.19 


_ 


_ 


14-141 


-1.07 


-1, 


.08 


-1, 


.04 






-1, 


.04 






-1.15 


_ 


_ 


11-2146 


-1.32 


-1 


.30 


-0 


.98 


-0, 


.75 


-1, 


.59 


-0, 


,68 


-1.20 


0.41 


-1.29 


11-1209 


-1.91 


-1, 


.92 


-1 


.74 


-1 


.70 


-1, 


.45 


-1, 


,60 


-1.95 


-1.23 


-0.22 


11-2322 


-1.53 


-1, 


.45 


-1 


.37 


-1 


.14 


-1 


.52 


-1, 


,07 


-1.12 


-0.48 


-0.54 


11-2037 


-1.75 


-1, 


.58 


-1, 


.25 


-1, 


.06 


-1, 


.78 


-0, 


.90 


-0.94 


0.05 


-1.25 


14-11 


-2.14 


-2 


.30 


-2 


.08 


-1, 


.27 


-2 


.44 


-1, 


.11 


-2.91 


0.89 


3.13 


14-125 


-1.43 


-1 


.33 


-1 


.32 






-1, 


.35 






-0.96 






11-2131 


-1.12 


-1 


.24 


-1 


.19 


-1 


.20 


-1 


.59 


-1, 


,04 


-1.70 


-1.12 


-1.53 


11-2031 


-1.31 


-1 


.17 


-0 


.90 


-0 


.73 


-1 


.46 


-0, 


,69 


-0.62 


0.20 


-1.77 


14-160 


-1.12 


-1, 


.00 


-0, 


.83 


-0, 


.65 


-1, 


.47 


-0, 


,60 


-0.54 


0.09 


-0.78 


14-183 


-1.52 


-1, 


.63 


-1, 


.53 


-1, 


.22 


-2, 


.08 


-1. 


,02 


-2.05 


-0.36 


-1.40 


CCDM-l-5734 


-3.04 


-3 


.06 


-2 


.80 


-2 


.15 


-2 


.74 


-1, 


.85 


-3.15 


-0.21 


-1.14 


21374275 


-0.81 


-0 


.71 


-0 


.74 


-0 


.65 


-1, 


.28 


-0, 


.59 


-0.35 


-0.52 


-0.70 


12-2519 


-1.70 


-1 


.55 


-1 


.42 


-1 


.32 


-1 


.80 


-1, 


,20 


-0.98 


-0.84 


-1.66 


12-1091 


-1.76 


-1 


.51 


-1 


.21 


-0 


.96 


-1 


.67 


-0, 


,89 


-0.53 


0.21 


-0.90 


13-1238 


-1.51 


-1 


.26 


-1 


.12 


-0 


.98 


-1 


.26 


-0, 


,87 


-0.25 


-0.40 


-0.21 


82-272 


-1.32 


-1, 


.33 


-1, 


.17 


-0, 


.92 


-1, 


.73 


-0, 


,71 


-1.36 


-0.05 


-1.31 


13-669 


-1.60 


-1, 


.39 


-1, 


.25 


-1, 


.16 


-1, 


.40 


-1, 


,11 


-0.59 


-0.68 


-2.10 


21380979 


-1.74 


-1, 


.54 


-1, 


.27 


-1, 


.09 


-1, 


.25 


-1, 


.03 


-0.75 


-0.15 


-0.39 


54-1781 


-1.90 


-1, 


.79 


-1, 


.54 


-1, 


.20 


-1, 


.68 


-1, 


,18 


-1.38 


0.04 


-1.67 


13-1877 


-1.42 


-1 


.25 


-1 


.01 


-0 


.80 


-1 


.31 


-0, 


,75 


-0.61 


0.14 


-1.80 


13-277 


-0.81 


-0 


.69 


-0 


.56 


-0, 


.48 


-1, 


.07 


-0, 


,51 


-0.20 


-0.04 


-1.24 


13-236 


-1.46 


-1 


.39 


-1 


.35 


-1, 


.38 


-1, 


.57 


-1, 


,31 


-1.10 


-1.37 


-1.61 


12-2113 


-1.42 


-1, 


.15 


-0, 


.96 


-0 


.77 


-0 


.94 


-0, 


,76 


-0.08 


0.05 


-2.09 


13-157 


-1.17 


-1, 


.22 


-1, 


.07 


-0 


.92 


-1, 


.63 


-0, 


,91 


-1.42 


-0.28 


-1.33 


13-52 


-2.22 


-1 


.97 


-1 


.59 


-1 


.43 


-2 


.15 


-1, 


,36 


-0.96 


-0.28 


-2.77 


91-155 


-1.12 


-0 


.92 


-0 


.69 


-0 


.50 


-0, 


.68 


-0, 


,51 


-0.15 


0.39 


-0.73 


21384350 


-1.76 


-1 


.67 


-1 


.51 


-1 


.77 


-1 


.64 


-1, 


,45 


-1.33 


-1.97 


-0.66 


54-1547 


-1.53 


-1, 


.48 


-1, 


.22 


-0, 


.88 


-1, 


.85 


-0, 


,74 


-1.32 


0.38 


-0.85 


91-506 


-1.73 


-1, 


.71 


-1, 


.46 


-1, 


.24 


-1, 


.81 


-1. 


,00 


-1.61 


-0.24 


-0.84 


12-1617 


-1.67 


-1, 


.66 


-1, 


.33 


-1, 


.09 


-2 


.00 


-0, 


.97 


-1.62 


0.11 


-1.68 


13-1250 


-1.72 


-1, 


.51 


-1, 


.14 


-0 


.76 


-1, 


.74 


-0, 


,76 


-0.66 


0.82 


-1.74 


21392541 


-2.12 


-1 


.83 


-1, 


.52 


-1 


.27 


-1 


.87 


-1, 


,23 


-0.69 


-0.10 


-1.13 


21392570 


-1.90 


-1 


.73 


-1 


.28 


-0 


.99 


-2 


.08 


-0, 


,85 


-1.08 


0.59 


-0.59 


21393104 


-1.63 


-1 


.58 


-1 


.31 


-1 


.07 


-2 


.09 


-0, 


,92 


-1.38 


0.01 


-1.15 


24-515 


-2.09 


-1, 


.98 


-1, 


.65 


-1, 


.42 


-1, 


.99 


-0, 


,97 


-1.55 


-0.24 


0.36 


21-998 


-1.77 


-1, 


.48 


-1, 


.42 






-1, 


.50 


-1, 


,19 


-0.33 




1.36 


21-33 


-1.69 


-1, 


.58 


-1, 


.41 


-1, 


.26 


-1, 


.41 


-1, 


.33 


-1.16 


-0.57 


-1.67 


92-393 


-2.00 


-1, 


.90 


-1, 


.72 


-1, 


.53 


-2, 


.06 


-1, 


,17 


-1.50 


-0.75 


0.96 


21395813 


-1.05 


-1 


.23 


-1 


.29 


-1 


.33 


-1 


.70 


-1, 


,28 


-1.93 


-1.56 


-1.77 


93-720 


-0.97 


-0, 


.86 


-0, 


.71 


-0, 


.67 






-0, 


,67 


-0.41 


-0.25 


-1.54 
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Table 11 — Continued 



Name 


a(H-5.8) 


a(H-8) 


a(H-12) 


a(H-17) 


q;(3.6-5.8) 


q;(H-24) 


a(5.8-8) 


a(8-17) 


q;(24-33) 


24-1796 


-1.97 


-1.90 


-1.48 


-1.13 


-2.02 


-0.75 


-1.63 


0.48 


-0.24 


21-2006 


-1.38 


-1.18 


-1.38 


-1.04 


-1.47 


-0.93 


-0.39 


-0.76 


-0.58 


KUN-196 


-1.72 


-1.89 


-1.99 


-1.44 


-2.86 


-1.29 


-2.58 


-0.48 


-0.29 


21402192 


-2.05 


-1.86 


-1.70 


-1.57 


-1.70 


-1.32 


-1.11 


-0.96 


-1.35 


23-570 


-1.53 


-1.48 


-1.38 


-1.28 


-1.48 


-1.20 


-1.29 


-0.84 


-1.23 


23-162 


-1.04 


-0.96 


-0.89 


-0.81 


-1.31 


-0.75 


-0.65 


-0.50 


-1.82 


23-798 


-1.35 


-1.31 


-1.05 


-0.72 




-0.72 


-1.15 


0.53 


-0.98 


01-1152 


-1.87 


-1.50 


-1.31 


-1.17 


-1.04 


-1.27 


-0.07 


-0.46 


-1.36 


01-580 


-1.50 


-1.59 


-1.38 


-1.13 


-1.81 


-1.05 


-1.94 


-0.15 


-0.90 


DG-481 


-2.02 


-1.96 


-1.80 


-1.88 


-1.82 


-1.92 


-1.72 


-1.72 


-2.75 



Note. — SED slopes (estimated as [log(\iFxl) — log(X2Fx2)]/(logXi — logX2)) used to analyze the disk struc- 
ture. The slopes are calculated using preferentially the IRS data (see text), and the IRAC and/or MIPS data for 
wavelengths not included in the IRS range and/or in case the IRS data presented a problem for particular objects. 
H photometry is taken from 2MASS. 



Table 12. Properties of the prototype disk models 



Model 


Mdisfc (M©) 


Rdisk (AU) 


Grain min-max, in/out (/im) 




Tin (K) 


Other / Comments 


Normal (1) 


0.015 


450 


0.1-10000 


0.16 


1500 


^cvys — 0.15 -AXJ 


Massive (2) 


0.052 


450 


0.1/0.1-100 


0.2 


1500 


rcomp— 3.0 AU 


Normal/Depleted (3A) 


0.00013 


100 


0.1-10000 


0.13 


1500 


Fits general structure 


Normal/Depleted (3B) 


0.00013 


100 


0.1/0.1-10000 


0.13 


1500 


rcomp=0.1 AU, improves near-IR excess 


Settled/Depleted (4) 


0.005 


200 


6-10000 


0.02 


1500 


Uniform grain composition 


Kink (5A) 


0.02 


150 


0.1/0.1-100 


0.08 


1500 


rcomp=2.0 AU; Good for silicate 


Kink (5B) 


0.02 


300 


10-10000/0.1-100 


0.08 


1500 


rcomp=5.0 AU; Good for general structure 


Old CTTS disk (6) 


0.0006 


200 


0.1/1-10000 


0.09 


1500 


rcomp— 2.0 AU 


Classical TD (7) 


0.0009 


200 


1-10000/0.1 


0.1 


1500 


rcomp— 5.0 AU 


Classical TD, sil (8A) 


0.0004 


100 


10-10000 


0.35 


250 


Good for general shape, fails to produce silicate feature 


Classical TD, sil (SB) 


0.0004 


200 


0.1/10-10000 


0.05 


150 


rcomp=6.0 AU; Probably needs a bit more dust in the hole 


Settled TD, no sil (9) 


0.003 


200 


20-10000/0.1 


0.04 


1500 


rcomp~15.0 AU 


Settled TD, sil (10) 


0.005 


200 


0.1/10-10000 


0.03 


1500 


rcomp=0.01 AU 



Note. — Properties of the RADMC models for the disk prototypes: disk mass and radius, maximum and minimum grain size, vertical pressure scale height at the outer disk radius 
(Rrdisk) in units of disk radius, inner rim temperature, and other properties (radius at which the dust composition changes [rcomp], radius below which all grains are crystalline 
[rcri/s]). The global dependency of the vertical pressure scale height goes as Hp/RocR^/'''. The power law exponent of the grain distribution is always assumed to be coUisional 
(a=-3.5), unless in the cases where a single grain population is assumed for part of the disk. The Schmidt number regulating the mixing of the two populations after the radius at 
which grain properties change (jcrys or rcomv) is always taken to be 2/3 (see Pavlyuchenkov & DuUemond 2007). 
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Table 13. Correlation table: Spearman rank coefficients 



Parameters 


R 


P 


R Range 


Comments 


Peak vs. Sp.Type 


-0.37 


0.006 


-0.61 0.06 


Excluding A,B stars 


Size vs. Peak 


-0.57 


9x10-5 


-0.77 0.26 




Peak vs. Age 


0.21 


0.14 


-0.11 — 0.50 


Excluding G stars or earlier 


M vs. Peak 


0.14 


0.46 


-0.31 — 0.54 




M vs. Size 


-0.10 


0.67 


-0.59 — 0.44 




q(H-5.8) vs. Peak 


0.27 


0.04 


-0.04 — 0.53 




a(H-8) vs. Peak 


0.28 


0.03 


-0.02 — 0.54 




a(H-12) vs. Peak 


0.48 


1x10""' 


0.20 — 0.68 




a(H-17) vs. Peak 


0.53 


6x10-5 


0.25 — 0.72 




a(5.8-8) vs. Peak 


0.33 


0.01 


0.03 — 0.58 




a(8-17) vs. Peak 


0.38 


0.005 


0.07 — 0.63 




q(H-24) vs. Peak 


0.37 


0.006 


0.06 — 0.61 




a(24-33) vs. Peak 


-0.41 


0.002 


-0.64 0.11 




a(H-5.8) vs. F/E 


0.33 


0.04 


-0.04 — 0.61 




a(H-8) vs. F/E 


0.27 


0.09 


-0.10 — 0.57 




a(H-12) vs. F/E 


0.13 


0.40 


-0.24 — 0.47 




a(H-17) vs. F/E 


0.16 


0.33 


-0.23 — 0.51 




0(5.8-8) vs. F/E 


-0.01 


0.93 


-0.37 — 0.35 




a(8-17) vs. F/E 


-0.12 


0.45 


-0.48 — 0.26 




a(H-24) vs. F/E 


0.17 


0.31 


-0.22 — 0.50 




a(24-33) vs. F/E 


-0.11 


0.49 


-0.46 — 0.27 





Note. — Correlation coefficients (R) and false-alarm probability (p) for different stellar, 
disk, and silicate properties, including height of the normalized peak, spectral types, 
size of the amorphous grains, age, forsterite to enstatite ratio (F/E), and SED slope at 
different wavelengths (a). 



